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ABSTRACT 

We use observed rotation velocity-luminosity (VL) and size-luminosity (RL) relations to single out a specific 
scenario for disk galaxy formation in the ACDM cosmology. Our model involves four independent log-normal 
random variables: dark-halo concentration c, disk spin Agai, disk mass fraction m„.^i, and stellar mass-to-light 
ratio T/ . A simultaneous match of the VL and RL zero points with adiabatic contraction requires low-c 
halos, but this model has V2.2 ^ l-SVvir (where V2,2 and Vvir are the circular velocity at 2.2 disk scale lengths 
and the virial radius, respectively) which will be unable to match the luminosity function (LF). Similarly 
models without adiabatic contraction but standard c also predict high values of V2.2/Vvii- Models in which 
disk formation induces an expansion rather than the commonly assumed contraction of the dark-matter halos 
have V2.2 ~ 1 •2yvir which allows a simultaneous fit of the LF. This may result from non-spherical, clumpy gas 
accretion, where dynamical friction transfers energy from the gas to the dark matter This model requires low 
Agai and mgai values, contrary to naive expectations. However, the low Agai is consistent with the notion that 
disk galaxies predominantly survive in halos with a quiet merger history, while a low mgai is also indicated by 
galaxy-galaxy lensing. The smaller than expected scatter in the RL relation, and the lack of correlation between 
the residuals of the VL and RL relations, respectively, imply that the scatter in Agai and in c need to be smaller 
than predicted for ACDM halos, again consistent with the idea that disk galaxies preferentially reside in halos 
with a quiet merger history. 

Subject headings: galaxies: formation — galaxies: fundamental parameters — galaxies: spiral — galaxies: 
structure 



1. INTRODUCTION 

In the standard, cold dark matter (CDM) based model for 
disk formation, set out by Fall & Efstathiou (1980), disks 
form out of gas that slowly cools out of a hot gaseous halo, 
associated with the dark matter potential well, while main- 
taining its specific angular momentum. During this process 
the dark matter halo contracts to conserve its adiabatic invari- 
ants (Blumenthal et al. 1986). Because of the centrifugal bar- 
rier the gas settles in a rotationally supported disk whose size 
is proportional to both the size and angular momentum of the 
dark matter halo (Mo, Mao & White, 1998; hereafter MMW). 
Consequently, the structure and dynamics of disk galaxies are 
expected to be strongly correlated with the properties of their 
dark matter halos. In particular, the correlations between the 
observable, structural parameters of disk galaxies, rotation ve- 
locity, V, size, R, and luminosity, L, are expected to be a re- 
flection of the virial properties of dark matter halos, which 

1/3 
scale as Vvii- oc /?vir « M^,.^; . Slight deviations from these scal- 

ings are expected from the fact that more massive halos are, on 

average, less concentrated. Any further deviations must either 

reflect some aspects of the baryonic physics related to galaxy 

formation, or signal a failure in the standard picture outlined 



above. In what follows we refer to the relations between the 
global disk parameters as the VL, RL, and RV relations. 

In the past, the VL relation, also known as the Tully-Fisher 
relation (Tully & Fisher 1977), has received much attention 
as a distance indicator owing to the relatively small observed 
scatter. Although numerous studies have addressed the origin 
of the VL relation, no consensus has been reached. In partic- 
ular, it is currently still under debate whether the origin of the 
VL relation is mainly governed by initial cosmological con- 
ditions (e.g., Eisenstein & Loeb 1996; Avila-Reese, Firmani 
& Hernandez 1998), or by the detailed processes governing 
star formation (Silk 1997; Heavens & Jiminez 1999) and/or 
feedback (e.g., Kauffmann, White & Guiderdoni 1993; Cole 
et al. 1994; EHzondo et al. 1999; Natarajan 1999). In the most 
recent models (van den Bosch 2000; 2002; Navarro & Stein- 
metz 2000; Firmani & Avila-Reese 2000, hereafter FAOO) it 
is typically understood that both initial conditions and bary- 
onic physics related to star formation and feedback must play 
an important role. Reproducing the VL zero point has also 
been a long standing problem for CDM based galaxy forma- 
tion models. In particular no model has been able to simul- 
taneously match the luminosity function and VL zero point 
using standard ACDM parameters (Cole et al. 2000; Benson 
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et al. 2003; Yang, Mo, & van den Bosch 2003). This prob- 
lem can be traced to the high values of V/Vvii expected for 
ACDM halos once the effects of the baryons, such as adia- 
batic contraction (Blumenthal et al. 1986), are taken into ac- 
count. All solutions to this problem require a change to either 
the standard cosmological model or the standard picture of 
galaxy formation. 

Another potential problem for disk formation models is the 
formation of large enough disks. The relative fragility of disk 
galaxies and the strongly ordered motion of their stars and gas 
is generally interpreted as evidence for a relatively smooth 
formation history without violent merger processes. The sub- 
sample of dark matter halos without recent major mergers is 
known to have systematically low spin parameters (D'Onghia 
& Burkert 2004). However, standard models for the formation 
of extended rotating disks seem to require high spin parame- 
ters in order to reproduce the zero point of the RL relation. 

In addition to the slopes and zero points of the VL and RL 
relations, additional constraints come from the scatter in these 
two relations. In particular, Courteau & Rix (1999; hereafter 
CR99) have shown that the residuals of the VL relation, at 
fixed L, are virtually uncorrected with the residuals of the 
RL relation, at fixed L (see also McGaugh 2005). This is an- 
other way of expressing the fact that the VL relation is inde- 
pendent of surface brightness (see Courteau et al. 2006 for 
details), and implies that the size of a disk at given L has no 
relevance to its rotation velocity. This is a puzzling result to 
explain, as one would naively expect that a more concentrated 
disk also results in a higher rotation velocity. CR99, there- 
fore, interpreted the weak residual correlation as indicating 
that, on average, high surface brightness (HSB) disks are sub- 
maximal, so that the disk only contributes mildly to the ob- 
served rotation velocity. If confirmed this puts constraints on 
the stellar mass- to-light ratio, T, and thus on the stellar initial 
mass function (IMF). FAOO, on the other hand, claim that the 
weak residual correlation owes to the surface density depen- 
dence of star formation, such that at a given baryonic mass, 
lower surface density galaxies have lower stellar masses, and 
hence lower luminosities, which compensate for the some- 
what lower rotation velocities. However, in their model both 
the highest and lowest surface brightness galaxies lie above 
the mean VL relation, contrary to observations, and their mod- 
els do not reproduce the zero point of the VL relation or the 
amount of scatter in the VL and RL relations. 

Although various studies have attempted to explain the ori- 
gin of the VL or the RL relation, the true challenge lies in find- 
ing a self-consistent model of disk formation that can match 
both relations as well as the galaxy luminosity function si- 
multaneously. Finding such a model is a non-trivial task, as 
all current models fail to do so. In this paper we examine the 
parameter space of such models within the standard ACDM 
cosmology. We simultaneously match the slopes, zero points 
and residuals of the VL and RL relations. We also investi- 
gate what each of these models predict for the mean V/Vvii- 
In order to be able to reproduce the observed abundances of 
disk galaxies, this ratio needs to be relatively low ^ 1.2. We 
show that this restriction severely limits the allowable param- 
eter space, favoring a model with low spin parameter Agai, low 
galaxy mass fraction mgai, and with halo expansion rather than 
halo contraction. 

The key ingredients of our model are as follows: 

1. Disk galaxies form in spherical dark-matter halos 
whose properties are drawn from N-body simulations 



of the standard ACDM cosmology. In particular, their 
density profiles have an NFW form with a concentra- 
tion parameter c that declines systematically with mass. 

2. The baryons form an exponential disk in centrifugal 
equilibrium, which is specified by a mass fraction mg^i 
and a spin parameter Agai. 

3. The galaxy mass fraction is treated as a free parameter 
with a mean njgai ex M^" . Positive values of a^ are 
expected from feedback effects. 

4. A bulge is included based on a self-regulating mecha- 
nism that ensures disk stability. 

5. The interaction between the baryons and the dark mat- 
ter halo is modeled with a generalized adiabatic con- 
traction model which also allows for halo expansion. 

6. Stars form in the disk once above a threshold surface 
density. 

7. The I-band stellar mass-to-light ratio T/ increases with 
luminosity as constrained by observations. 

8. The model parameters c, Agai, wigai and T/ are assumed 
to be independent log-normal random variables. 

The data and models are introduced in ^ and ^ respec- 
tively. In 33.4l we outline the conversion between stellar mass 
and luminosity; In ^3.51 we detail the computation of our 
model scaling relations. In ^we discuss how to construct 
models that match the slopes, zero points, scatter and residual 
correlation. In ^we explore the model parameter space, and 
advocate a revised model for disk formation. We summarize 
our results in ^ 

2. THE DATA 

We compare our models with the large data set of ~ 1600 
local disk galaxies compiled by Courteau et al. (2006). For 
each of these galaxies a rotation velocity V , a luminosity L, 
a disk scale length R, and a central surface brightness of the 
disk /zo, are available. The data set is compiled from three 
independent samples: Mathewson, Ford, & Buchhorn 1992 
(hereafter, MAT); Courteau et al. 2000 (hereafter, Shellflow); 
and Dale et al. 1999 (hereafter, SCII). 

All galaxy luminosities are measured in the /-band, and cor- 
respond to total luminosities (i.e., disk plus bulge). In addi- 
tion, for a significant sub-sample we also have optical colors 
avail able, which we use to estimate stellar mass-to-light ratios 
(see ^3.4l below). Disk scale lengths and surface brightnesses 
are determined from the /-band photometry. In what follows, 
whenever we refer to surface brightness we mean the central 
surface brightness, /xqj, of the exponential disk fitted to the 
data. 

Due to the complications involved with interpreting Hi line 
widths, we only use rotation velocities derived from resolved 
Ha rotation curves. This reduces the full sample by ~ 300 
galaxies. For the MAT and Shellflow samples the rotation 
curves are fitted with a parametric function (Courteau 1997) 
which is then evaluated at 2.2 disk scale lengths. For the 
SCII sample the velocities are measured at the optical radius 
(equivalent to 3.2 scale lengths for an exponential disk). 
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Fig. 1 . — Observed I-band VLR scaling relations using data from Courteau et al. 2006. Bi-weighted orthogonal least squares fits are given by the solid black 
lines, with 2-sigma deviations given by the dashed lines. The open black circles with en'or bars show the mean and 2-sigma scatter of the VL and RL relations 
binned at 0.3 dex intervals in L/. The colors and point types correspond to extrapolated disk central surface brightness, figj as indicated in the top right panel. 



2.1. Corrections 

In order to homogenize these data samples as much as pos- 
sible, we have derived the inclination and cosmological cor- 
rections to the velocities, luminosities and scale lengths in a 
uniform way, as described below. 

The observed rotation velocities, Vobs, are corrected for in- 
clination and cosmological broadening using 



V: 



Kbs 



(l+z)sin/ 
The inclination, /, is computed using 



sm; : 



-ib/a)^ 



1 



-?o 



(1) 



(2) 



where b/a is the minor-to-major axis ratio, and ^o is the in- 
trinsic thickness of the disk. We assume ^o = 0.2, and set 
i = 90° if b/a <qo. 

The apparent magnitudes, mj, are corrected for both inter- 
nal extinction, Aint, and external (i.e.. Galactic) extinction. 



Aext- A small k-correction, Ak, is also applied such that 

mi = m/,obs -Ai„t-Aext-Ak. (3) 

Internal extinctions are computed using the line-width (W = 
2V) dependent relation from TuUy et al. (1998): 

Ain, = -f!(W)log(a/b) = [0.92+ 1.63(logH'- 2.5)] log(a/b) 

(4) 
with a/b the major- to-minor axis ratio. The external (Galac- 
tic) extinction is computed using the dust maps of Schlegel, 
Finkbeiner, & Davis (1998), while the k-corrections are com- 
puted using the line width dependent formalism of Willick 
etal. (1997). 
The absolute magnitudes, M/, are computed using 



M/ = m/-51ogDL-25; Dl = 



VcMB 

100/j 



(l+z), 



(5) 



with VcMB the systemic velocity of the galaxy in the reference 
frame at rest with the cosmic microwave background (Kogut 
et al. 1993). The /-band luminosities are computed from Mj 
using an absolute magnitude for the Sun of M/ q = 4.19. 
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Fig. 2. — Residual correlations from the observed I-band VLR relations in Fig.[T] The black lines show weighted least-squares fits. The RL residuals show a 
clear correlation with surface brightness, while the VL residuals show none. The residuals of the VL and RL relation are only weakly correlated. The dotted fine 
has a slope of -0.5 expected for a pure exponential disk (CR99). 



Disk scale lengths are corrected for inclination using 

^i = -Ri,ob.s/[l+0.41og(fl/fo)], (6) 

(Giovanelli et al. 1994), and converted into kilo parsecs using 
the angular diameter distance Da = Vcmb/[100/i (1 +z)]. 

Finally, central surface brightnesses are corrected for incli- 
nation. Galactic extinction, and cosmological dimming {per 
unit frequency interval) using 

Mo,/ = Mo,/,obs + 0.51og(a/^)-Aext-2.51og(l+z)l (7) 

The factor of 0.5 in front of the \og{a/b) term is empirically 
determined by demanding that the residuals of the relation be- 
tween central surface brightness and rotation velocity has no 
inclination dependence'. Following Giovanelli et al. (1997) 
we assume an uncertainty of 15% in 7/, a/b, Agxt, and A^, and 
propagate the errors. Not including distance uncertainties, the 
average errors on the observables are ainv — 0.08, uinL — 0.1, 
and (TinR ~ 0.14. 

' For a disk of zero thickness one expects the factor to be between 0, for 
an optically thick disk, and 2.5, for an optically thin disk. 



2.2. Scaling Relations 

We now use the data described above to investigate the VL 
and RL relations of disk galaxies. Despite the fact that our 
sample consists of subsamples that use slightly different tech- 
niques or definitions for the scale lengths and rotation veloci- 
ties, we find that, with the uniform inclination and extinction 
corrections described above, each subsample yields VL and 
RL relations that are consistent with each other (see Courteau 
et al. 2006 for details). We therefore combine the three sub- 
samples to a single sample of ~ 1300 galaxies. 

The resulting VLR relations are shown in FiglT] The solid 
lines show the mean relations, which have been determined 
using bi-weighted orthogonal least-squares fits. The dashed 
lines show the 2a scatter. The open circles with error bars 
show the data mean and la scatter in separate luminosity bins 
with a width of 0.3 dex. These show that the VL and RL re- 
lations are well fitted with a single power-law over the range 
considered here. The RV relation has the largest scatter, and 
is thus the least well determined. For consistency we use the 
mean VL and RL relations to determine the mean RV relation. 
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Fig. 3. — Comparison of our I-band VLi? scaling relations with those from the literature. The long-dashed lines show the VL relation from Giovanelli et al. 1997 
(G97), and the dotted lines show the RY relation from Courteau (1996, 1997). The red lines are derived from Pizagno et al. (2005). 



The mean relations thus derived are: 

log ^^ = 0.296 log -J^^ 

^ [kms-i] ^ [10"'-3/272Lq] 

log— i^^^ — = 0.3221og ^' 



[/z^^kpc] 



[1010-3/2-2Lq] 



+ 2.162, (8) 
+ 0.455. (9) 



log 



R, 



[/lyokpc] 



1.086 log 



y 



[kms-i] 



1.894. 



(10) 



where Ri is the de-projected disk scale length in the photo- 
metric /-band, and /ivo = //o/(70kms"' Mpc~'). 

The color coding in FiglT]denotes the central surface bright- 
ness of the disk in the /-band, /xq./. Note that the scatter in 
the YL relation is not correlated with /iq./. This is also ap- 
parent from Fig. |2] which shows the correlations between the 
residuals of the YL and RL relations at constant L (AlogV 
and Alog/?, respectively) and surface brightness, /io/. Again 
the symbols are color coded according to /io,/. As is evident 
from the upper left panel, the scatter in the RL relation is by 



definition^ dominated by scatter in surface brightness. The 
residuals of the YL relation, however, show no significant cor- 
relation with surface brightness at all (lower left panel). In 
addition, the YL and RL residuals are only very weakly corre- 
lated (lower right panel): the slope of the residual correlation 
is 7 = d[ A log y(L)] /d[ A log/?(L)] = -0.08 ±0.03, in agree- 
ment with CR99. 

FigElshows a comparison of the YLR relations derived here 
with previous studies. The black lines and point types are the 
same as in Fig.[T] The long-dashed and dotted lines show the 
YL relation from Giovanelli et al. (1997), and the RM relation 
that fits the data of Courteau (1996, 1997). Note that these 
YL and RL relations, which are in excellent agreement with 
our mean scaling relations, were used by MMW and FAOO as 
model constraints. The red lines in Fig. [3] show the YLR re- 
lations of Pizagno et al. (2005). The latter study is based on 
a sample of 81 disk dominated (B/D < 0.11) galaxies, with 

^ A family of pure exponential disks with constant T; has, at a constant 
Z,/, AlogR oc Q.ifiQj, which is different from the observed slope of 0.13. 
This owes to the relatively small range of luminosities sampled, relative to 
the amount of scatter in the RL relation. 
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luminosities in the Sloan i-band, disk scale lengths measured 
with bulge-to-disk decompositions, and velocities measured 
at 2.2 disk scale lengths from Ha rotation curves. We convert 
from Li to Lj using logL/ = logL, +O.4(/-/)-O.4(/-/)0 ~ 
0.036, assuming M,.© = 4.56, and (/-/) ~ 0.46 (Courteau 
et al. 2006). Note that the slope of the VL relation of Pizagno 
et al. (2005) is somewhat steeper than ours, while their RL 
relation is offset from ours towards larger disk scale lengths. 
These differences may be due to a combination of their rela- 
tively small sample size (81 galaxies), their bulge-to-disk ratio 
selection criteria, and different inclination corrections. These 
differences do not significant impact on our main conclusions, 
as we discuss in AppendixlAl 

3. DISK FORMATION MODELS 

In the standard picture of disk formation (e.g.. Fall & Efs- 
tathiou 1980), disks form inside virialized dark matter halos 
through the cooling of the baryonic material. Our models are 
based on this standard picture, and closely follow MMW, but 
with some additional ingredients. In particular, we include 
a model for the formation of bulges, we use a prescription 
for star formation which separates the disk into gaseous and 
stellar components, we use an improved and generalized de- 
scription for adiabatic contraction, and we use empirical rela- 
tions to convert our models to observable quantities for direct 
comparison with the data described above. All these model 
ingredients are discussed in more detail below. 

Although we refer to our models as 'disk formation mod- 
els' they are completely 'static' (i.e., the model is not actu- 
ally evolved). For a given specific angular momentum of the 
baryonic material out of which the disk forms and a given 
potential due to the dark matter, the structural properties of 
the resulting disk are computed assuming (i) that disks are 
exponential, and (ii) that the specific angular momentum of 
the baryonic material is conserved. Alternatively, one could 
in principle consider more 'dynamic' models, that follow the 
actual formation of the disk galaxies starting at high redshifts 
(e.g., Firmani & Avila-Reese 2000; Avila-Reese & Firmani 
2000; van den Bosch 2001, 2002). However, as long as the 
formation of the disk is sufficiently quiescent, the final struc- 
ture of the disk should be independent of its actual formation 
history: the structural properties of the final disk are basically 
just governed by the principle of dynamical, centrifugal equi- 
librium. In a more dynamic approach, one can actually model 
the star formation history of the disk, which is not possible 
with our static model. However, the star formation history 
mainly governs the final mass-to-light ratio of the stars, which 
we set using empirical relations. The models described here 
should thus be applicable independent of the detailed forma- 
tion history as long as the two assumptions mentioned above 
are satisfied. 

Unless stated otherwise we adopt a ACDM cosmology with 
n^ = 0.3, ilb = 0.044, r^A = 0.7, h = //o/(100kms-i Mpc"') = 
0.7 and with a scale-invariant initial power spectrum with a 
normalization as = 0.9. The baryonic mass fraction of this 
cosmology /bar — 0.15. 

3.1. Disk Formation 

We model dark matter halos as spheres with a NFW density 
distribution 

P^"^= . / V?I / ^2 (11) 

(r/rs)(l + r/r,y 

(Navarro, Frenk & White 1997) where r^ is a characteristic 
radius at which the logarithmic slope of the density distribu- 



tion dlnp/dlnr = -2, and ps = p('"s)- The overall shape of the 
density distribution can be characterized by the so-called con- 
centration parameter c = /?vir/'"s- Here Ryii is the virial radius, 
which is defined as the radius inside of which the average halo 
density is Ayir times the critical density for closure. For the 
ACDM cosmology adopted here Ayir ~ 100 (Bryan & Nor- 
man 1998). In addition to the virial radius /?vir we also define 
the virial velocity Vvii as the circular velocity at the virial ra- 
dius, i.e., y^^j. = GM^iy/Ryiy, with G the gravitational constant. 
The total angular momentum of a halo, 7vii, is commonly 
expressed in terms of the dimensionless spin parameter: 



A = 



GMf 



V 2^virVvir 



Here E is the halo's energy, and fc measures the deviation 
of E from that of a singular isothermal sphere with the same 
mass, and is given by 



c 1 - 1/(1 +c)2-21n(l +c)/(l +c) 
2 [c/(l+c)-ln(l+c)]^ 



(13) 



(see MMW). 

We assume that the galaxy that forms consists of a bulge 
and a disk. Our algorithm to ascribe a bulge-to-disk ratio is 
discussed at the end of this section. We define the baryonic 
mass of the total galaxy as 

Mgal = Md+Mb = mgalMvir (14) 

with Md and M\, the masses of the disk and bulge, respectively, 
and < Wgai < /bar- We define the bulge-to-disk mass ratio as 
9 = Mb/Md so that 

1 



Md = 



-m„alMvi 



Mb: 



i+e""^"' ™ 
e 

y-j-^mgalMvir. 



(15) 



(16) 



In addition, we write that the total angular momentum of the 
baryons out of which the disk plus bulge form is 



■'gal — JgiiWvi 



(17) 



As shown by van den Bosch et al. (2002), the specific angular 
momentum distribution of the total baryonic mass (including 
those baryons in the halo that do not partake in the formation 
of the disk plus bulge) is virtually identical to that of the dark 
matter Therefore, one also expects that < jgai ^ /bar- 

As discussed below, we assume that the bulge forms out of 
disk instabilities. Let Jh = jhJvW indicate the original angular 
momentum of the baryonic material out of which the bulge 
forms. We assume, however, that the bulge formation process 
transfers this angular momentum to the disk plus the halo, so 
that the final angular momentum of the bulge is zero. Indeed, 
such an angular momentum transfer is observed in numerical 
simulations (e.g., Hohl 1971; Debattista et al. 2006). If we 
define / as the fraction of J\, that is transferred to the disk, we 
obtain that the final angular momentum of the disk is equal to 

^d=[;gai-(l-/t);b]/vir (18) 

The specific angular momentum of the final disk is therefore 



A = (i+e)[i-/uf^)^ 

where we have introduced the parameter 



/i 



lost '■ 



:(l-/t)( — 

^^gal 



(19) 



(20) 
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which expresses the fraction of the total angular momentum 
of the material out of which the bulge plus disk forms that has 
been lost to the halo. For modeling purposes it is more useful 
to define the parameter 



A = /lo 



1 + 9 

e 



■d-ft) 



Jb/Mb 



(21) 



^7gal/Mgai^ 

which expresses the ratio of the specific angular momentum 
that has been lost to the halo due to bulge formation to the 
total specific angular momentum of the material out of which 
the disk plus bulge have formed. For example, /x = 1 means 
that bulge formation does not change the specific angular mo- 
mentum of the disk. The extreme /x = occurs when the disk 
loses mass but not angular momentum during bulge forma- 
tion. Note that unlike /lost, the parameter /x can in principle 
be larger than unity. In practice, however, the bulge is likely 
to form out of material with relatively low specific angular 
momentum (e.g., Norman, Sellwood & Hassan 1996; van den 
Bosch et al. 2002). Furthermore, the fraction of angular mo- 
mentum that is lost to the halo is expected to be fairly small 
(e.g., Weinberg 1985; Debattista & Sellwood 2000; Valen- 
zuela & Klypin 2003; O'Neill & Dubinski 2003), so that fi is 
expected to be close to unity. 

For a disk with a surface density T,{R) and a circular veloc- 
ity V(R) the total angular momentum is 

Jd = 2TT / j:(R)RV(R)RdR. (22) 



Throughout we assume that the disk that forms has an expo- 
nential surface density distribution T,{R) = ^oexp{-R/Rd), so 
that 

yd = 2Md/?dKir/v (23) 



with 



J l-Rvir/Ra 



h = \l e"« 



.„..2VXM«d) 



dM 



(24) 



(cf., MMW). If we combine equations ([T2]l, (O, and ( l23T l 
we obtain the following expression for the scale length of the 
disk: 

7?d= ^4^ [l + (l-/x)e] AgaiT?™ (25) 



/vV2A 

where we have defined Agai as the effective spin parameter of 
the material out of which the bulge plus disk form: 



A, 



gal : 



7gal 
'"gal 



A 



(26) 



From equation ( |25] | it is evident that bulge formation impacts 
the final scale length of the disk. Typically, for /x > 1 the disk 
size will decrease, while /x < 1 causes an increase in /?d- Note 
that the transition does not occur exactly at /x = 1, because the 
presence of a bulge component modifies /y. Numerical simu- 
lations indicate that bulge formation causes an increase in disk 
scale lengths (Debattista et al. 2006), suggesting that /x < 1. 
Shen et al. (2003) considered a fairly similar model but with 
the assumption that the material that forms the bulge has the 
same specific angular momentum as the disk. Their favored 
model is equivalent to our model with /x = 0.5. However, 
given the arguments above, we expect /x to be smaller than 
this. In what follows we adopt a fiducial value of /x = 0.25, 
although none of our results are very sensitive to this particu- 
lar choice. 

Note that since the computation of /y requires knowledge 
of Ri, this set of equations has to be solved iteratively (see 



MMW for details). When computing the total circular veloc- 
ity V{R) we use that 



Here 



V\R) = V^{R) + V^{R) + Vly,{R) (27) 



yliR) = ^2/ [Io(y)Ko(y)-h(y)Ki(y)] 

"d 



(28) 



with y = R/(2Rd), and /„ and K„ are modified Bessel functions 
(Freeman 1970). The circular velocities of the bulge, Vh, and 
the dark matter halo, Vdm, are computed assuming spherical 
symmetry, whereby the mass distribution of the dark matter 
halo is adjusted for adiabatic contraction (see ^3.21 below). 
The density distribution of the bulge is assumed to follow a 
Hernquist (1990) profile 



Pb- 



Mb 



rb 



2tt r(r+rb) 



(29) 



In projection, this is similar to a de Vaucouleurs profile (i.e., a 
Sersic profile with Sersic index n = 4), with a half light radius, 
Ref{= 1.8152rb (Hernquist 1990). Throughout we adopt the 
relation between R^g and M\, from Shen et al. (2003): 



log/?eff = 



-5.54 + 0.561ogMb (logMb > 10.3) ,„, 
-1.21+0.141ogMb (logMb<10.3) ^^^' 



Although bulges of late-type disk galaxies are better described 
by exponential profiles (e.g. Courteau, de Jong & Broeils 
1996), what matters most for our purposes is the total bulge 
mass; its distribution is only of secondary importance. 

The computation of the actual bulge-to-disk ratio for each 
model galaxy is based on the fact that self-gravitating disks 
are unstable against global instabilities (e.g., bar formation). 
We follow the approach of van den Bosch (1998, 2000) and 
Avila-Reese & Firmani (2000) and assume that an unstable 
disk transforms part of its disk material into a bulge compo- 
nent in a self-regulating fashion, such that the final disk is 
marginally stable. Bars, which are considered the transitions 
objects in this scenario, are thus expected to be fairly com- 
mon, as observed. 

We define P(R) = Vd(R)/V(R), and consider the disk to be 
stable as long as 



Pn 



■■ max P(R)<Pcnt 

0<R<Rvir 



(31) 



(Christodoulou, Shlosman & Tohline 1995). The actual value 
of /3crit depends on the gas mass fraction of the disk, but falls 
roughly in the range 0.52 < f3„h Js 0.70. For a given value 
of /3crit, we use an iterative technique to find the bulge-to-disk 
ratio for which /9max = /^crit- 

3.2. Adiabatic Contraction 

When baryons cool and concentrate in the center of a dark 
matter halo they deepen and modify the shape of the gravi- 
tational potential. If this process is slow with respect to the 
(local) dynamical time of the halo, the halo will contract to 
conserve its adiabatic invariants. For the idealized case of 
a spherical halo in which all dark matter particles move on 
circular orbits, the adiabatic invariants reduce to the specific 
angular momentum, rV(r). If, in addition, the distribution of 
the baryons has spherical symmetry, this reduces further to 
rM(r), with M(r) the enclosed mass within radius r (Blumen- 
thal et al. 1986; hereafter BFFP). 

In realistic dark matter halos, however, the particles typ- 
ically move on highly eccentric orbits (Ghigna et al. 1998; 
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Fig. 4. — Ratio of observed to virial circular velocities, V2.2/Vvir> as a 
function of the concentration parameter, for various forms of adiabatic con- 
traction. All models have Ajai = 0.048 and nigai = 0.05. The vertical dotted 
lines show the mean c according to model of Bullock et al. (2001a) for ba- 
les with virial mass, Mvi, = 10'^,10'^, & 10"/i"'Mq. The horizontal dotted 
lines show VS j/l^vir ratios of 1 to 1.8 at intervals of 0.2. The relation using 
the standard (Blumenthal et al. 1986) adiabatic invariant, rM{r), is given by 
the blue sohd hne. The modified adiabatic invariant, rM{f), as proposed by 
Gnedin et al. 2004 results in only a ~ 0.02 dex (~ 5%) reduction in V2.2 (red 
short dashed line). Most of this reduction is taken back if we use specific an- 
gular momentum as the adiabatic invariant, and take into account the disk ge- 
ometry when computing Kirc (dot-dashed red line). For all of these adiabatic 
contraction models, standard concentration parameters yield Va.^/K-ir — 1-6- 
The black long-dashed lines show models where we have artificially modi- 
fied the contraction factor (see equation [34l . Starting from the standard model 
(v = 1), the rM(r) relation is approximately reproduced with v = 0.8, v = Q 
results in no adiabatic contraction of the halo, while v =—\ gives halo ex- 
pansion. The green dot-long dashed line shows Vmax/K-ir. where Vniax is the 
maximum halo circular velocity (without adiabatic contraction). 



van den Bosch et al. 1999). Taking this into account reduces 
the effect of adiabatic contraction (Wilson 2003). Gnedin 
et al. (2004) give a modified adiabatic invariant, rM(r), where 
r ~ 0.85(r//?vir)"'^ is the orbit averaged radius. Using this adi- 
abatic invariant results in somewhat less contraction of the 
halo than in the standard BFFP formalism. An additional 
problem is that disks are not spherical. Therefore, in principle 
one should use rV(r) as the adiabatic invariant, rather than 
rM{r). It is well known that the circular velocity curve of 
a thin exponential disk rises less rapidly but reaches a higher 
peak velocity than a sphere with the same enclosed mass (e.g., 
Binney & Tremaine 1987). Therefore, using rVir) as an adi- 
abatic invariant, rather than rM(r), results in a stronger con- 
traction of the dark matter halo at 2.2 disk scale lengths. 

In addition to these somewhat subtle problems for the stan- 
dard BFFP formalism, one may also question whether adia- 
batic contraction really occurs during disk formation. If disks 
are not built by smooth, relatively slow, spherical infall, the 
adiabatic contraction could, in principle, be counter-balanced 
by a variety of processes. These may even go as far as to cause 
an actual expansion of the dark matter distribution. One such 
process is dynamical friction (e.g. El-Zant, Shlosman & Hoff- 
man 2001; Ma & Boylan-Kolchin 2004; Mo & Mao 2004; 
Tonini, Lapi & Salucci 2006), which could be significant if 
disks are built by relatively big clumps. Indeed, the physics 
of gas cooling suggests that disks may have formed out of 
clumpy, cold streams rather than from a smooth cooling flow 
(Birnboim & Dekel 2003; Mailer & BuHock 2004; Dekel & 



Birnboim 2006; Kaufmann et al. 2006; Keres et al. 2005). 
We defer a detailed study of the effects of such cold infall on 
the contraction of the halo to a future paper. Here we con- 
sider a simple modification of the BFFP adiabatic contraction 
formalism that allows us, with a single tunable parameter, to 
consider reduced contraction, no contraction, or even expan- 
sion. 

Our method starts from the BFFP formalism, according to 
which a dark matter particle initially (i.e., before the forma- 
tion of the disk plus bulge) at radius ri settles at a radius rf, 
where 

rfMf(rf) = nM,(n). (32) 

Here Mi{r) and Mf{r) are the initial and final mass distribu- 
tions. If we assume that initially the baryons have the same 
(normalized) density distribution as the dark matter, then 
Mi(r) is simply given by the initial halo profile (e.g. NFW). 
For the final mass distribution we have that 



Mf (rf) = Md(rf ) + Mb(rf) + Mqm j('-f) 

= Md(?-f) +Mb(rf) + ( 1 - mgai)Mi(n) 



(33) 



Here Md(r) is the mass of the exponential disk enclosed within 
spherical shells of radius r, Mb(r) is the similarly enclosed 
mass of the Hernquist bulge, and the second equality follows 
from the assumption that adiabatic contraction occurs with- 
out shell crossing. Equations ( |32] | and ( l33T l can be solved it- 
eratively for the contraction factor r(ri) = rf/ri, which then 
allows for a computation of the mass distribution of the con- 
tracted dark matter halo. 

Our modification consists of simply defining the actual re- 
lation between rf and r; as 



rf=r''ri 



(34) 



with v a free parameter: v = I yields the standard BFFP con- 
traction, 1^ = corresponds to no adiabatic contraction, and 
ly < Q models an expansion of the dark matter halo. As a 
specific example of an expansion model, v = -\ results in an 
expansion factor that is equal to the contraction factor in the 
BFFP model. 

Fig. m shows the impact of various adiabatic contraction 
formalisms on the ratio of the total circular velocity at 2.2 
disk scale lengths, V2.2, to the virial velocity, Vvir. We con- 
sider models with an effective spin parameter Agai = 0.048 
and a baryonic mass fraction m^^\ = 0.05 (no bulge forma- 
tion is considered here). We apply the various contrac- 
tion formalisms described above and compute the resulting 
V2.2/Vvir as a function of the halo concentration parameter 
c. The three vertical dashed lines (from left to right) in- 
dicate the mean halo concentration expected for halos of 
mass Mvir = 10'^ 10'-, &IO" /i-'M© in the model of Bullock 
etal. (2001a). 

Note that with the standard BFFP adiabatic contraction 
(v = 1), and a mean halo concentration parameter for a halo 
of mass Mvii- = 10'^/i"'Mo, we expect that Va.a/Kir — 1-7, 
while V2.2/Kir — 1-45 without adiabatic contraction {v = 0). 
In the extreme case where the halo is adiabatically expanded 
{v = -1) the ratio is further lowered to V2.2/Vvir — 1-3. As an 
illustration we also plot Vmax/Kir, where Vmax is the maximum 
circular velocity of a NFW halo, which is related to the halo 
concentration as 



Vr. 



Vvi 



' 0.465 



(35) 



ln(l-l-c)-c/(l+c)' 
Note that for the typical concentration of galaxy sized halos 
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The short dashed (red) line in Fig.|4]shows the results for the 
adiabatic invariant rM{f) proposed by Gnedin et al. (2004). 
As mentioned above, this adiabatic invariant results in a some- 
what smaller overall contraction than the standard BFFP for- 
malism. Note that we can accurately model the Gnedin 
et al. formalism by simply setting v = 0.8. Not only does this 
reproduce V2.2 but also the shape of the circular velocity pro- 
file. Finally, the dot-short dashed (red) line shows the results 
for the adiabatic invariant rV{r), with F the orbit averaged ra- 
dius of Gnedin et al. (2004). This adiabatic invariant accounts 
for both the eccentricity of the orbits of the dark matter par- 
ticles and for the (non-spherical) geometry of the disk. Note 
that the resulting Vaa/Kh is virtually indistinguishable from 
what one obtains with the standard BFFP formalism: taking 
account of the non-sphericity of the disk completely cancels 
the impact of non-circular orbits. 

3.3. Star Formation 

Disks are made up of stars and cold gas. In the high- 
est surface brightness galaxies the gas fraction (defined as 
the ratio of cold gas mass to total disk mass) is small, typ- 
ically ~10%, so assuming a pure stellar disk is reasonable. 
However, the gas fraction increases with decreasing surface 
brightness and luminosity (McGaugh & de Blok 1997; Kan- 
nappan 2004) such that low surface brightness galaxies have 
comparable amounts of mass in cold gas and stars. The data 
used here, however, only provides measurements of the stel- 
lar disks, while the models describe the distribution of total 
baryonic matter (cold gas plus stars). To allow for a proper 
comparison between models and data we need a prescription 
for computing the ratio between stars and cold gas as function 
of radius in the galaxy. 

Kennicutt (1989) has shown that star formation is strongly 
suppressed below a critical surface density, which can be 
modeled by a simple Toomre stability criterion: 

E„,(/;)=^Iil^ (36) 

3.36 gG 

(Toomre 1964). Here K{R)=^/2'^{\ + i^|^)^- is the epicy- 
cle frequency, cTgas is the velocity dispersion of the gas, and 
Q is the Toomre stability parameter. Throughout we adopt 
CTgas = 6kms"' and Q=\.5; with this choice of parameters the 
radius at which star formation is observed to truncate coin- 
cides with the radius where Sgas = Sciit (Kennicutt 1989). 
We compute the total stellar mass as 

M, =Mb + 27r/ [•>:{R)-i:^ni{Ry\RAR (37) 

JO 

where 7?sf is the star formation truncation radius, defined by 
S(/?sf) = Sci-it(/?sF)- Note that we thus assume that the bulge 
is made of stars entirely, and that all the disk material with a 
surface density above the critical density has been converted 
into stars. Although this is clearly an over-simplification of 
the complicated physics associated with star formation, there 
are several reasons why this is probably reasonable. Firstly, it 
is a clear improvement over the assumption that disks consist 
of stars only, such as in the standard MMW model. Secondly, 
as shown by van den Bosch (2000, 2001), this simple model 
(i) matches the gas mass fractions of disk galaxies as function 
of their surface brightness and luminosity, and (ii) naturally 
leads to gas disks that are more extended than stellar disks. 
Finally, more detailed, 'dynamic' models for disk formation, 
which model the actual star formation rate using realistic, em- 
pirically motivated, prescriptions, show that the typical star 



formation time scale is short compared to the time scale on 
which the disk accretes new gas. Consequently, these mod- 
els indeed predict that the gas disk has been depleted by star 
formation down to Sgas = Sent (van den Bosch 2001). As we 
will see later, the inclusion of a star formation threshold den- 
sity proves a crucial ingredient in solving two problems of the 
standard MMW model. 

3.4. Conversion from Mass to Light 

In order to compare our models with observations, we 
need to convert both the stellar masses and the stellar disk 
scale lengths into the observed /-band luminosities and scale 
lengths, respectively. 

The conversion from mass to light is conventionally done 
via the mass-to-light ratio, T = M,/L. Bell & de Jong (2001) 
showed that T can be estimated from optical colors. These re- 
lations have been updated by Bell et al. (2003a) and Portinari, 
Sommer-Larsen & Tantalo (2004). The main uncertainty in 
this method is the normalization, reflecting the unknown stel- 
lar IMF. Additional uncertainties arise from details related to 
the amount of dust extinction and the star formation histo- 
ries. Upper limits on the normalization can be obtained from 
maximal disk fits to observed rotation curves, as shown in 
Bell & de Jong (2001), although these may still suffer from 
distance uncertainties. More accurate estimates from rotation 
curve mass modeling are hindered by well known degenera- 
cies (e.g., van Albada et al. 1985; van den Bosch et al. 2000; 
Dutton et al. 2005). 

For the sub-sample of our galaxies with optical colors, we 
apply the relations in Bell et al. (2003a), to obtain estimates 
of T/. These are based on a 'diet'-Salpeter IMF, introduced 
by (Bell & de Jong 2001). We have verified that, assuming 
the same IMF, the relations of Portinari et al. (2004) give very 
similar results that agree to within 0.05 dex. In Fig|5]we plot 
the resulting T/ as a function of the /-band luminosity. The 
solid red line shows the bi-weighted orthogonal least squares 
fit to these points, and is given by 



log 



T, 



[Mo/Lq] 



= 0.172 + 0.144 log 



[W^%ILq] 



+ A 



IMF, 



(38) 

The intrinsic scatter in this relation is estimated to be ~ 0.1 
dex (Bell et al. 2003a, Kauffmann et al. 2003a). We model 
this scatter by drawing T/ for each model disk galaxy from 
a log-normal distribution with a mean give by eq. ( [38] l and 
with a scatter cinx, which we treat as a free parameter. The 
other free parameter, Aimf, absorbs our lack of knowledge 
about the IMF. It is equal to zero for the 'diet'-Salpeter IMF 
assumed here, while Aimf — +0.15 for a Salpeter (1955) IMF, 
and Aimf — -0.15 for a top-heavy IMF like that of Kennicutt 
(1983). See Portinari et al. (2004) for details. 

As a consistency check, we apply the Bell et al. (2003a) 
relations to the mean (g - r) color-luminosity relation of Piza- 
gno et al. (2005). The resulting T/(L/) is indicated as the 
long dashed line in Fig. |5] and agrees with our mean rela- 
tion to within 0.05 dex. For comparison, we also show the 
T/ adopted by previous studies. The short-dashed line shows 
the relation adopted by FAOO; the relatively weak slope of 
their T/(L/) relation was derived from a comparison of the 
slopes of /-band and //-band Tully-Fisher relations, while 
their normalization is based on the somewhat ad hoc assump- 
tion that the disk contributes 70 percent to the circular veloc- 
ity at 2.2 disk scale lengths. Finally, McGaugh et al. (2000) 
and MMW both adopted a constant mass-to-light ratio (in- 
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Fig. 5. — I-band stellar mass-to-light ratio, T/, versus I-band luminosity. 
The points show T; estimated from optical colors using the relations in Bell 
et al. (2003a). These data form a sub-sample of the Courteau et al. (2006) 
data set: MAT (green). Shell (blue), and SCII (red). A typical Icr error bar 
is given in the lower right corner The solid red line gives the bi-weighted 
fit to these data. As a consistency check the long dashed line gives the rela- 
tion computed using the mean (g — r) color -luminosity relation from Pizagno 
et al. (2005) and the T,-(g - r) color relation from Bell et al. (2003a). We 
transform /-band luminosities into /-band luminosities with L/ = L, -0.036 
(Courteau et al. 2006). This relation is in excellent agreement with ours. For 
comparison we give the T; = 1.7 estimate of McGaugh et al. (2000; dot- 
dashed line), T/(L/) relation from FAOO (short-dashed line), and the adopted 
T/ = l.lh (dotted line) and the used T; = 1.19/i (dot-long dashed line) by 
MMW. 



dependent of luminosity). The value adopted by McGaugh 
et al. (2000), T/ = 1.7(M/L/)q, is based on a stellar popula- 
tion synthesis model with a Salpeter IMF, while MMW based 
their value, T/ = \.lh(M /Li)q on the sub-maximal disk ar- 
guments of Bottema (1993). Note, however, that there is 
an error in MMW, and that the actual value they used is 
T/ = 1.19Ii(M/Lj)q (H.J. Mo 2004, private communication). 
Note that this value is small compared to the typical mass-to- 
light ratios shown here. As we will demonstrate below, this 
has important implications for the MMW results. 

3.4.1. Disk Scale Lengths 

In order to compute the /-band scale lengths of our model 
disks, for direct comparison with the data discussed in ^ we 
proceed as follows. 

For each galaxy we first compute the stellar surface density 
profile of the disk, i.e., E, = T,(R)-J:„JR) with 0<R< Rsf, 
which we fit with an exponential disk over the radial range 

0.l5RsF<R<0MRsF. 

In theoretical models of disk galaxy formation, disks form 
inside out. This results in color gradients, with progressively 
larger scale lengths when going from stellar mass to /T-band 
light (assumed to closely trace the underlying stellar mass) to 
B-band light. Observations show a similar trend with Rr/Rh = 
1.17 (MacArthur, Courteau, & Holtzman 2003), and Ri/Rk = 
1.13 (de Jong 1996), where subscripts denote the photometric 
band. To account for these color gradients we convert our 
stellar disk scale lengths to /-band scale lengths using log/?/ = 
log/?* -I- 0.05, with /?* the scale length of the exponential fit 
to the model stellar disk. Note that this assumes the scale 
length in stellar mass, /?», is equal to the scale length in K- 



band light, Rk, and thus is probably an underestimate of the 
true correction. 

3.5. Sampling Strategy 

In order to pursue a meaningful comparison of our model 
VL and RL relations with the observations, we construct sam- 
ples of random realizations of model galaxies. Each model 
galaxy is specified by four parameters; the virial mass of the 
halo, Mvii, the halo concentration parameter, c, the effective 
spin parameter of the galaxy, Agai, and the baryonic mass frac- 
tion of the galaxy, mgai. 

Previous studies (e.g., MMW; Somerville & Primack 1999; 
Cole et al. 2000; van den Bosch 2000; Firmani & Avila-Reese 
2000; Croton et al. 2006) have assumed that Agai = A. How- 
ever, there are numerous reasons why the spin parameter of 
the disk may be different from that of its dark matter halo (see 
©. In this paper we assume that Agai follows a log-normal 
distribution: 



/'(Agal)dAgal=' 



1 



flnA^ 



:exp 



ln^(Agal/Agal) 



2-L 



dAgal 



gal 



(39) 



and we treat Agai and ^inA as free parameters. If Agai — A we 
expect that Agai — 0.042 and a\n\ ~ 0.5, which are the values 
obtained for dark matter haloes from cosmological numerical 
simulations (Bullock et al. 2001b). In what follows we will 
consider these as our fiducial values. 

The concentration parameter, c, is strongly correlated with 
the halos mass accretion history (Wechsler et al. 2002; Zhao 
et al. 2003; Li, Mo & van den Bosch 2005), and thus depends 
on both halo mass and cosmology. Bullock et al. (2001a) and 
Eke, Navarro, & Steinmetz (2001) present analytical models, 
calibrated against numerical simulations, for the computation 
of a mean concentration given a halo mass and cosmology. In 
what follows we use the model by Bullock et al. (2001a). At a 
given halo mass, the halo concentrations follow a log-normal 
distribution with a scatter crinc = 0.32 (Wechsler et al. 2002). 
However, as for the spin parameter, the mean and scatter are 
different for the subset of halos without a recent major merger. 
As shown by Wechsler et al. (2002), the scatter becomes 
smaller while the mean increases. To be able to account for 
this, we consider ainc a free parameter and compute c(Mvir) 
using the relation by Bullock et al. (2001a) but multiplied by 
a free parameter rj^ ■ 

Unlike for c and A, very little is known regarding the bary- 
onic mass fractions, mgai, of galaxies. We account for this 
limitation by modeling the mean as 



Wlgal(Mvir) = mgaLO 



M„i 



lOii-5/i-iM, 



(40) 



with mgai,() and a^ two free parameters. The value of a,n is re- 
lated to the relative efficiencies of cooling and feedback pro- 
cesses. Typically, cooling results in am < while feedback 
results in Ofm > (e.g., Dekel & Silk 1986; van den Bosch 
2002). Finally, to allow for scatter in the above relation, we 
assume that, at fixed halo mass, m^ai follows a log-normal 
distribution with scatter crin,,,, which we also consider a free 
parameter 

To construct VL and RL relations for our models we proceed 
as follows. We uniformly sample a range in logMvii. For each 
halo, we then draw values for c, Agai, and njgai using the log- 
normal distributions described above. We then iterate until a 
solution for the disk-to-bulge ratio is found, taking adiabatic 
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TABLE 1 

Overview of Model Parameters and their Fiducial 

Values 



Parameter 


symbol 


fiducial value 


Median effective spin parameter 


Aeal 


0.042 


Concentration parameter normalization 


^c 


1.0 


Adiabatic contraction parameter 


V 


1.0 


Toomre stability parameter 


Q 


1.5 


Bulge formation thi'eshold 


/Sen. 


1.0 


Bulge formation exchange parameter 


/x 


0.5 


Galaxy mass fraction normalization 


"Igal.O 


0.05 


Galaxy mass fraction slope 


am 


0.0 


Mass-to-light ratio normalization 


AlMF 


0.0 



contraction into account, and compute the resulting galaxy lu- 
minosity, Li, disk scale length Ri, (both in the /-band), as well 
as the circular velocity, V2.2, of the model galaxy at 2.2/?/. To 
mimic the observational errors we add to each of the three 
model observables logV22, logL/, and log/?/ a Gaussian de- 
viate with a dispersion that reflects the typical observational 
errors. The uniform sampling of halos in logMvir results in ap- 
proximately uniform sampling in logL/, while the luminosity 
sampling of the data is approximately log-normal with a low 
luminosity tail. To allow for a fair comparison with the data 
we need to reproduce the observed luminosity sampling. We 
do this with a Monte-Carlo technique which accepts model 
galaxies such that they reproduce the observed luminosity dis- 
tribution. Finally, when including bulge formation, we tune 
Pa-it SO as to roughly reproduce the observed distribution of 
bulge-to-disk ratios of our data sample. 

3.6. Overview of Model Parameters 

A list of the free model parameters with their fiducial values 
is given in Table 1. The fiducial value Pa-it =1-0 corresponds 
to a model without bulge formation. 

In addition to these model parameters, there are four param- 
eters that describe the amounts of scatter in the log-normal 
distributions of c, Agai, 'Wgai and T/, and which we also treat 
as free parameters in what follows. 

4. COMPARISON BETWEEN MODELS AND DATA 

In order to gain some insight into the origin of the slopes 
and zero-points of the VL and RL relations, we start by con- 
sidering a set of simplified models without bulge formation 
(/3crit = 1 -0), and without scatter {ux^c = crinA = o-inm = crinx = 0). 
More complete models, including both bulge formation and 
scatter will be presented in ^4.21 

4. 1 . Median Parameter Relations 

4.1.1. Slopes 

We start by reproducing the models of MMW. To this ex- 
tent we set am = 0, 2 = 00 (i.e., entire disk is made of stars), 
and T/ = 0.83 (i.e. 1.19 x 0.7). As can be seen in the left- 
hand panels of Fig. |6] this model reproduces the slope of the 
VL relation, but predicts a somewhat steeper slope for the RL 
relation than observed. The deviations of the slopes of both 
the VL and RL relations from the virial value of 1 /3 are due 
to the non-homology of the dark matter halos (i.e., the c(Mvii) 
relation). Our deviations from the pure virial equations are 
somewhat larger than in MMW. This owes to the fact that we 
use the halo concentration model of Bullock et al. (2001a), 
which predicts a much stronger mass dependence than the 
NFW model used by MMW. The Buflock et al. (2001a) model 
also predicts higher c than the NFW model. This accounts for 



the slight discrepancy between the VL zero point of our model 
and the data. 

Panels in the second column from the left show the same 
model, but with the more realistic T(L) relation of eq. ( |38] ). 
Note that the slopes of both the VL and RL relations are now 
significantly steeper than the data. In the third column of 
Fig.|6l we include the star formation threshold to separate the 
disk into stars and gas. Since lower luminosity galaxies have 
higher gas-to-stellar mass ratios this flattens the VL slope, 
bringing it back in agreement with the data-'. The star forma- 
tion threshold also results in smaller stellar disk scale lengths 
compared to the baryonic disk scale lengths. However, this 
reduction is approximately canceled out by the reduction in 
luminosity so that the RL slope and zero points are not sig- 
nificantly affected by the star formation threshold. Finally, in 
order to match the RL slope we tune «„, to match the slope of 
the RL relation. This requires a^ = 0.18, so that the galaxy 
mass fraction, mgai, increases with increasing halo mass, as 
expected, for example, for simple supernova feedback models 
(e.g., Dekel & Silk 1986; van den Bosch 2002; Dekel & Woo 
2003). As shown in the right-hand panels of Fig.|6]this model 
roughly matches the slopes of both the VL and RL relations, 
and will hereafter serve as our reference model. 

Although the introduction of a non-zero «„ causes the 
galaxy mass fractions, m„.ji, to systematically vary by a fac- 
tor ~ 3 over the luminosity range probed, this does not sig- 
nificantly impact the slope of the VL relation. This owes to 
the fact that variance in mgai scatters galaxies mainly along 
the VL relation: an increase in Wgai makes the galaxy more 
luminous but simultaneously increases its rotation velocity 
(see also Navarro & Steinmetz 2000). This important aspect, 
which we discuss in more detail in 34.2l below. implies that we 
can simply tune a^ and ;«gai.o to fit the slope and zero point of 
the RL relation, without (strongly) affecting the VL relation. 

4.1.2. Zero points 

Having explored how some of our model parameters impact 
the slopes of the VL and RL relations, we now turn to the zero 
points, which we define as the values of V and R at the mean 
luminosity of our data; logL/ = 10.3. For our reference model 
with am = 0.18 which matches the slopes of the VL and RL 
relations, the model galaxies are both too large (by ~ 35 per- 
cent) and rotate too fast (by ^ 40 percent). This corresponds 
to a 2(7 offset in terms of the observed scatter in the VL rela- 
tion, and thus an even more significant offset in terms of the 
actual uncertainty in the VL zero point. 

The inability of CDM based galaxy formation models to 
match the VL zero point is a generic problem that has been 
identified in numerical simulations (Elizondo et al. 1999; 
Navarro & Steinmetz 2000, Eke, Navarro & Steinmetz 2001; 
Sommer-Larsen et al. 2003), in semi-analytic models (e.g.. 
Mo & Mao 2000; van den Bosch 2000, Firmani & Avila- 
Reese 2000; Cole et al. 2000; Benson et al. 2003) and in halo 
occupation models (Yang, Mo, & van den Bosch 2003). In 
particular no model has been able to simultaneously match 
the luminosity function and the zero point of the VL rela- 
tion, using standard ACDM parameters. Cases that claim an 

' We emphasize that the luminosity dependence of the gas-to-stellar mass 
ratio is the key ingredient to this success, not the mechanism by which this 
is achieved. Therefore any argument against S^rit as a physical threshold for 
star formation is not relevant here. All that matters is that our separation of 
disk mass in stars and cold gas reproduces the gas-to-stellar mass ratios as 
a function of disk luminosity and surface brightness, as is the case for our 
models. 
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Fig. 6. — VL and RL relations using our fiducial mean c, Xg^i, and m^.^i parameters and different values of T/, Q and a,,,. The black solid and dashed lines show 
the mean and 2cr scatter of the observations. The red line is a fit to the model galaxies, whose slope is given in the top left comer of each box. The color of the 
points corresponds to /lo.; as in Fig.[T] The model on the far left is the MMW mod el with T/ = 0.83. Q = oo (a pure stellar disk) and Qm = 0. The model in the left 
middle only differs from this model in the T;, which we assume follows Eq.|38] This steepens the slopes of the VL and RL relations making them incompatible 
with the data. The right middle model includes the star formation threshold, which results in a gas fraction that increases with decreasing luminosity, and hence 
a shallower VL slope, in agreement with the data. The effect on the RL slope is not as strong because as the luminosities are decreased, so are the stellar disk 
scale lengths. Notice that the star formation threshold has introduced a small curvature to both the VL and RL relations. In order to match the RL slope we set 
Om = 0.18 (far right panel). Although m„.j[ now varies significantly with luminosity, the only effect on the VL relation is to increase the curvature slightly. While 
this model matches the slope of the VL and RL relations it fails to reproduce the zero points. 



agreement either assume that the observed rotation velocity, 
Vobs, is equal to Vvir (e.g., Somerville & Primack 1999), or to 
Vmax (e-g- Croton et al. 2006). In both cases, the effect of 
the baryons on the rotation curve is completely ignored. As 
shown in Fig. HI when adiabatic contraction of the dark mat- 
ter halo and the disk's own contribution to Vob,s are properly 
accounted for, one predicts that Vobs/Kir — 1-7 for a typical 
halo in a ACDM cosmology. Models that fail to take these 
effects into account can therefore not be used for a meaning- 
ful comparison with the data (see also Navarro & Steinmetz 
2000; van den Bosch 2000). For example, Croton et al. (2006) 
model Vobs as the maximum circular velocity of a NFW halo, 
Vmax- As shown by the dot-dashed curve (labelled Vmax/Vvir) 
in Fig.m this more or less corresponds to halo expansion with 
i/ = -1 .0. We will return to the implications of this in ^ 

In principle, there are a number of different ways in which 
one might envision solving the VL zero point problem: 

1. Lower stellar mass-to-light ratios. This trivially in- 
creases the luminosity of the model galaxies at a fixed 
rotation velocity. For our reference model we need to 
lower T by ~ 0.5 dex. However, the most that one can 
justify based on realistic IMFs is Aimf — -0.2 dex (cor- 
responding to a Kennicutt IMF), more than a factor of 
2 smaller than what is required. 

2. Lower halo concentrations. This results in a lower V2.2 
at a given L. Halo concentrations can be lowered by, for 
example, decreasing the cosmological parameters il,^ 
and/or erg . A reduction of the power spectrum on small 
scales also results in lower halo concentrations (Zent- 
ner & Bullock 2002). However, for our reference model 
to match the VL zero point we need rjc ~ 0.3 which is 



TABLE 2 
Model Parameters 



Model 


A,Mf 


V 


■r]c 


e 


-^gal 


'Wgal.O 


dm 


An. 


/x 


I 


-0.40 


0.8 


0.80 


1.5 


0.042 


0.06 


0.25 


0.66 


0.25 


II 


-0.20 


0.8 


0.50 


1.5 


0.042 


0.10 


0.25 


0.75 


0.25 


III 


-0.20 


0.0 


0.80 


1.5 


0.042 


0.09 


0.25 


0.80 


0.25 


IV 


-0.20 


-1.0 


1.00 


1.5 


0.023 


0.03 


0.30 


0.85 


0.25 



difficult to reconcile with current constraints on cosmo- 
logical parameters and the matter power spectrum. 



3. Modify adiabatic contraction. As discussed in ^3.2l a 
more realistic treatment of adiabatic contraction than 
the standard BFFP formalism is unlikely to have any 
significant impact on V22- However, if we simply turn 
off adiabatic contraction (i.e. v = 0), V2.2 will be low- 
ered by ~ 20%. An even stronger reduction can be ac- 
complished by considering expansion, modelled by set- 
ting i^ < in Eq. ( l34b . 

Based on these possibilities we construct three models that 
match the VL and RL zero points and slopes. The parameters 
of these models are listed in Table 2. In each of these models 
we tune the parameter mgai and a^ to fit the RL zero point 
and slope. Note that these models should only be considered 
specific examples of a more extended parameter space, which 
we describe in more detail in ^ These models mainly serve 
to highlight the various possible solutions to the VL zero point 
problem discussed above. 

In model I we consider elements from all three modifica- 
tions discussed above. We reduce the halo concentrations by 
20% (i.e. we set rjc = 0.8), which corresponds to changing 
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(n^,cr^) from (0.3,0.9) to (0.25,0.8), as advocated by van 
den Bosch, Mo & Yang (2003) and which is in agreement with 
the third year WMAP results (Spergel et al. 2006). In addi- 
tion, we use the adiabatic invariant of Gnedin et al. (2004) 
which corresponds to i^ = 0.8. Finally, we adjust the stel- 
lar mass-to-light ratios until we match the zero-points, which 
requires Aimf = -0.4. As discussed above, such a large re- 
duction in T/ implies an unrealistically top-heavy IMF. This 
model should therefore be considered as an illustration only. 
Because of the lower c and the reduced halo contraction, the 
resulting disk galaxies are larger than in our reference model. 
We counter-balance this by slightly increasing mgai.o from our 
fiducial value of 0.05 to 0.06. Alternatively, we could have 
matched the increase in disk scale lengths by decreasing Agai: 
such models are discussed in ^ 

In model II we restrict Ajmf to -0.2, which is the most we 
can accommodate with realistic IMFs, and we match the VL 
zero point by lowering c by 50%. As shown by Zentner & 
Bullock (2002), such a large reduction in halo concentrations 
can be reconciled with a power spectrum with a running spec- 
tral index that is still consistent with the WMAP data (Spergel 
et al. 2003). 

Finally, in model III we simply turn off adiabatic contrac- 
tion (ly = 0). This model is able to match the VL and RL zero 
points with the standard c(Mvii) for a ACDM cosmology with 
n„, = 0.25 and cg = 0.8 (i.e., 77c = 0.8) and with a realistic IMF 
(Aimf = -0.2). 

4.2. Models with Scatter 

Having identified models that can simultaneously match the 
slopes and zero points of the VL and RL relations, we now 
focus on the scatter in both relations. The observed scatter 
in the VL and RL relations at fixed Lj are aiav = 0.13 and 
cins = 0.32 respectively. The observational uncertainties are 
estimated to be <t\„v — 0.08, uinL — 0.10, and crins ~ 0.14 (see 
^2.1l i. Subtracting these in quadrature (with an error weighted 
scheme) from the observed scatter leaves an intrinsic scatter 
of a\nv = 0.12 and crins = 0.28. Note that our error weighted 
subtraction results in a slightly larger intrinsic scatter than ob- 
tained from a straight quadratic subtraction. 

The scatter in our models originates from four sources: T/, 
c, Agai, and mg.^\. Table 3 lists the expected amount of scatter 
in each of these four variables, and Fig|7] shows the affect of 
these sources of scatter on the VL and RL relations. Lower 
c, higher Agai, and lower mgai all result in larger scale lengths, 
and lower stellar masses and luminosities. Higher T/ result 
in lower luminosities, and higher effective (i.e. at a fixed lu- 
minosity) scale lengths. The fact that scatter in c and Agai 
also effects the luminosity is entirely due to the star forma- 
tion threshold and results in a significantly reduced VL scatter 
due to Agai. As already discussed in ^4.1.11 scatter in mgai 
moves galaxies approximately parallel to the VL relation, but 
only for intermediate values of the galaxy mass fraction: for 
large mgai scatter results in V oc L"^, while for sufficiently low 
values of mgai the corresponding slope approaches zero. 

Fig.|2]also yields useful insight on the effect of the various 
parameters on the zero points of the VL and RL relations. The 
most effective way to change the VL zero point is to change 
the average halo concentration, c, or the average stellar mass- 
to-light ratio, T/. The zero point of the RL relation, on the 
other hand, is most easily changed by a modification of the 
mean spin parameter, Agai, or the average galaxy mass frac- 
tion, mgai. 



TABLE 3 
Sources of intrinsic scatter 

""In T, ""In c 0"ln A ""In mj 



Maximum 0.32 
Predicted 0.23 
Adopted 0.23 



0.40 
0.32 
0.23 



0.25 0.55 

0.50 

0.25 0.00 



Note. — The first row gives the maximum scat- 
ter in T/, c, Agai, and mg;,i that is consistent with the 
VL and RL scatter. The second row gives the pre- 
dicted scatter, for c and Agai from cosmological sim- 
ulations, for T/ from observations. The third row 
gives the scatter adopted in models I-III in order to 
match the VL and RL scatter simultaneously while 
including bulge formation. 

4.2.1. Constraining the amount of scatter 

To determine the maximum amount of scatter allowed for 
each parameter we run models with log-normal scatter in one 
of the four parameters only. We increase the scatter until we 
reach the intrinsic scatter of either the VL or RL relation. The 
resulting limits on the amounts of scatter in c, Agai, Wgah T/, 
are very similar for each of the three models discussed above. 
In Table 3 we list the representative values thus obtained. As 
expected from Fig|7] the limits for T/ and c come from the 
VL relation, while those for Agai and mgai owe to the RL re- 
lation. Note that the expected amounts of scatter (also listed 
in Table 3) in T/ and c are smaller, and therefore consistent 
with, their respective maximum amounts. In the case of the 
spin parameter, however, the expected amount of scatter is a 
factor two larger than the maximum amount allowed by the 
intrinsic scatter in the RL relation. This can have a number of 
important implications, of which we consider the following 
three: 

1. Due to surface brightness selection effects our data is 
missing the lowest and/or highest surface brightness 
galaxies. Although this is certainly possible to some 
extent, gauging from the incompleteness corrections of 
de Jong & Lacey (2000) we conclude that these effects 
are likely to be small (Courteau et al. 2006). 

2. Disk galaxies acquire a narrower distribution of specific 
angular momentum than their host halos (i.e., the scat- 
ter in In Agai is smaller than that in In A). This can occur, 
for example, through the redistribution of angular mo- 
mentum associated with bulge formation. We discuss 
such a possibility in the next section. 

3. Disk galaxies form in a sub-set of halos with a smaller 
CTin A- Interestingly, this is expected if disk galaxies form 
mainly in those halos that have not experienced recent 
major mergers (D'Onghia & Burkert 2004). In addition 
to a smaller scatter in both A and c, this subset of halos 
also has a significantly smaller mean spin parameter, 
and possibly a higher mean c (Wechsler et al. 2002). 
We discuss such a model in ^ 

Having determined the maximum amounts of scatter al- 
lowed by the data, we now proceed to set their actual val- 
ues. Since scatter in mj does not contribute much scatter to 
the VL relation, and the scatter in the RL relation is already 
over-budgeted because of the spin parameter, we set a\nm = 
in what follows. In addition, we set the scatter in the stellar 
mass-to-light ratios to the expected amount, i.e., crinx = 0.23. 
If we also set a\nc to its predicted value of 0.32, then the scatter 
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Fig. 7. — Contribution of scatter in T;, c, \„.^\, and m„.^\ to the I-band VL and RL scaling relations for our reference model with Myij = 3 X lO" /r' Mq 
and AiMF = ~0.5 (to match the zero points). Models with/without the star formation threshold are given by circles/triangles. The solid symbol shows the mean 
parameter model, the open symbols show the models with ± 1 & 2o- scatter in the parameter. The arrow in the top left of each panel shows the direction in which 
the parameter increases. The observed mean and 2a scatter are given by solid and dashed lines, respectively. See text for further details. 



in the VL relation is already larger than observed, even with- 
out any scatter in Agai. We therefore adopt cr\„c = 0.23, which 
is the value predicted for halos that have not experienced any 
recent major mergers (Wechsler et al. 2002). Finally, we tune 
the scatter in Agai until we match the observed scatter in the 
RL relation, which yields a\n\ = 0.25 (when including bulge 
formation, see below). Throughout this paper we treat c, T/, 
mgai and Agai as independent variables. In reality, the scatter in 
some of these parameters may be correlated. In Appendix IBJ 
we discuss possible correlations, and how they impact on our 
results. 

4.2.2. Bulge Formation 

Thus far we have only considered models without bulge for- 
mation (i.e., with f3cm = 1 in Eq. dSTI )). We now introduce 
bulge formation to the three models described above. We set 
/x = 0.25, and tune f3a-n so that the sample of model galaxies 
has the same mean bulge-to-disk ratio as the data (~ 0.15). 
None of our results are sensitive to the exact value of /x cho- 
sen. 

Since mgai increases and c decreases with halo mass, more 
massive halos typically host galaxies with a larger /3max- Also, 
for a given halo mass, the disk's contribution to the circu- 
lar velocity increases with decreasing spin parameter Agai. 
Consequently, it is predominantly the luminous, high surface 
brightness galaxies residing in massive halos that will form a 
bulge in our models. 

Bulge formation only has a small effect on the VL relation, 
but does cause a modest change of the RL relation. For our 
fiducial value of /x, and in fact for all /x ^ 1, bulge formation 
causes an increases of the specific angular momentum of the 
disk, and thus of its scale length. Since bulge formation pref- 
erentially affects the galaxies in the lower right part of the RL 
plane (i.e. the highest surface brightness galaxies), this will 
actually reduce the scatter in the RL relation at the bright end 
(see also Shen et al. 2003). At the faint end no bulges form. 



so that this also causes a small increase in the slope of the 
RL relation, which we counter by increasing a^. For realis- 
tic bulge-to-disk ratios the overall scatter in the RL relation is 
reduced at most by ~ 15%. Although this also changes the 
maximum amounts of scatter allowed by the data by a sim- 
ilar fraction, the maximum scatter in Agai allowed by the RL 
scatter is never larger than 0.3, which is still much smaller 
than the predicted amount. Thus bulge formation is not able 
to reconcile the expected scatter in A with that of the observed 
scatter in the RL relation. 

The VL and RL relations for models I-III with bulge forma- 
tion, scatter in c, Agai, T/, and observational errors are shown 
in Fig.|8] Overall, all three models provide a reasonable match 
to the slopes, zero points and scatter of the observed relations. 
A detailed comparison with Fig. [T] however, reveals that all 
three models, but especially Model 11, have a slight problem 
at the bright end where they predict rotation velocities that are 
somewhat too high. This is related to the fact that the contri- 
bution of the baryons to V2.2 increases with luminosity (see 
34.3.2b . Consequently, the VL relation deviates from a pure 
power-law. As we will see in 34.3l below. this causes a weak, 
but significant correlation between surface brightness and the 
VL residual. 

4.3. Residual correlations 

Models I, II and III all fit the slopes, zero points and the 
amounts of scatter of the VL and RL relations. An impor- 
tant question is whether we can discriminate between these 
three models, or whether there are genuine degeneracies in 
the model parameter space. As shown in Fig |2] there is ad- 
ditional information in the residual correlations. In particu- 
lar, the residuals in the VL and RL relations are virtually un- 
corrected, so that the VL residuals are not significantly cor- 
related with surface brightness (CR99). We now investigate 
how models I, II and III fare in matching these aspects of the 
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Fig. 8. — I-band VL and RL scaling relations for our three models with bulge formation and scatter in c, Ag^i, T; and observational errors. The parameters of 
these models are given in Table 2. The black solid and dashed lines show the mean and 2cr scatter of the data, respectively. The solid red lines give the mean 
relation of the model galaxies from unweighted least-squares fits of V on L and R on L. The colors and point types correspond to surface brightness as in Fig[T] 
Each model sample consists of ~ 1300 galaxies, sampled to reproduce the observed distribution of luminosities. 



data. 

The upper panels of Fig. |9] show the residuals of the RL re- 
lations of models I-III plotted against central surface bright- 
ness. All three models are virtually indistinguishable, and 
accurately reproduce the data (upper-left panel of Fig. |2]i. 
This, however, owes simply to the luminosity sampling of the 
data, which we reproduce using a Monte-Carlo technique (see 
^2.21) . The lower panels of Fig. |9] however, plot the VL resid- 
uals against central surface brightness. This time, the three 
models differ slightly from each other, and significantly from 
the data. All models predict a small but non-negligible cor- 
relation between surface brightness and the VL residual, con- 
trary to the data (lower-left panel of Fig.|2]i. This owes mainly 
to the non power-law character of the model VL relations, as 
discussed above. 

The upper panels of Fig.fTOlshow the VL residuals as func- 
tion of the RL residuals. Except for Model 11, which predicts 
a correlation slope 7 = d[ A log y(L)] /d[ A log/?(L)] which is 
somewhat too steep, these residual correlations are in good 
agreement with the data (lower-right panel of Fig. |2|i. 

As discussed in |T] the absence of a pronounced correlation 
between the VL and RL residuals is contrary to naive expec- 
tations. In order to understand which aspects of the model 
cause this success, we start by considering the model bary- 
onic VM and RM relations (i.e. velocity at 2.2 baryonic scale 
lengths vs. baryonic mass and baryonic scale length vs. bary- 
onic mass). The lower panels of Fig.fTOlshow that the residu- 
als of these baryonic scaling relations are strongly correlated, 
for each of the three models. This owes to the fact that a 



smaller Agai implies both a reduction of the disk scale length 
and an increase in Vi.i- Thus, the naive expectation that the 
residuals should be correlated at least holds for the baryonic 
scaling relations. We find this to be a very generic predic- 
tion of our models, even though the exact value of the slope 
of the baryonic residual correlation may vary somewhat from 
model to model. Although measuring baryonic masses and 
baryonic sizes depends on stellar mass-to-light ratios, which 
are difficult to constrain on a galaxy by galaxy basis, these 
can be determined in a statistical sense for a suitably chosen 
sample. Thus, in principle the baryonic scaling relations are 
observable, and our predictions regarding their residual corre- 
lations provide a useful test for our models. We caution, how- 
ever, that it is crucial that the rotation velocity is measured at 
2.2 disk scale lengths. For example, when V is measured at 
5 disk scale lengths, our models predict a significantly shal- 
lower residual correlation. 

The middle panels of Fig.fTOlshow the residual correlations 
from the model stellar mass VM and RM relations (i.e. ve- 
locity at 2.2 stellar scale lengths vs. stellar mass and stellar 
scale length vs. stellar mass). Note that the strength of these 
correlations is significantly reduced with respect to the bary- 
onic case (lower panels). This is caused by the star formation 
threshold, which correlates the stellar mass at a given bary- 
onic mass to the spin parameter: a lower Agai not only results 
in a smaller scale length and a larger V2.2, it also results in 
a relatively larger stellar mass, therewith reducing the resid- 
ual correlation. This, together with the constraint of match- 
ing the slope of the VL relation, is another argument in favor 
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Fig . 9 . — Residuals of the model /-band VL and RL relations versus surface brightness, for models in Fig.[8] The observed relations are given by the dashed line. 
All three models reproduce the surface brightness dependence of the RL relation, but none of the models are able to reproduce the observed surface brightness 
independence of the VL relation. 



of including a star formation threshold (or something equiva- 
lent) in our models. In fact, FAOO argued that a star formation 
threshold density is the main explanation for the weak residual 
correlation observed. Although we agree that it substantially 
reduces the residual correlations, it does not automatically re- 
sult in uncorrelated residuals, as is evident from the middle 
panels of Fig. [TOJ The weaker correlation between the VL and 
RL residuals in the upper panels compared with the middle 
panels in Fig. [TOJis achieved through scatter in T/ and, to a 
lesser extent, to our modelling of the observational errors in 
V, L and R. 

4.3.1. Impact of Scatter 

To gauge how the various sources of scatter impact the 
residual correlations, FigfTTI shows AlogV against AlogT? 
for models that have only scatter in one of the four param- 
eters: T/, c, Agai, or mgai. The spin parameter, Agai, and the 
galaxy mass fraction, mgai, only cause significant scatter in 
the RL relation. Consequently, they result in residual corre- 
lations with I7I < 0.1, in agreement with the data, which has 
7 = -0.08 ± 0.03 (see ^2.2b . Scatter in c and T/, however, 
predict strongly correlated residuals with 7 ~ -0.6 and 7 ~ 1, 
respectively. 

Thus, achieving VL and RL relations with uncorrelated 
residuals requires a subtle balance between the amounts of 
scatter in the halo concentrations, c, and the stellar mass- 
to-light ratios, T/. To emphasize this point, FiglT2] shows 
the residual correlations for models with different amounts 
of scatter in c and T/. The scatter in Agai and the scatter due 
to observational errors are kept constant. A model with only 
scatter in Agai and observational errors has 7 = -0.06. Includ- 



ing the maximum amount of scatter in c (middle panel) results 
in a strong negative correlation, 7 = -0.19. In contrast, includ- 
ing the maximum scatter in T/, with no scatter in c, results in 
a weak, positive correlation with 7 ~ 0.06 (right panel). In 
practice scatter in both c and T/ are expected to contribute to 
the VL scatter, so there is not much freedom to tune the resid- 
ual correlations. In order to reproduce the total VL scatter we 
adopted crinc = crinx = 0.23. As shown in Figs. |9]&[l0] this 
amount of scatter results in a weak residual correlation and a 
small but significant surface brightness dependence to the VL 
scatter. Removing these dependences entirely would require a 
larger scatter in T/ and a lower scatter in c. The former is fea- 
sible, but the latter would require an even more biased subset 
of halos than just those without recent major mergers. 

4.3.2. Dark and baryonic mass fractions 

The upper panels of Fig. [13] show the ratio Vgai/Vi.a as func- 
tion of surface brightness for all three models. Here Vgai and 
V2.2 are the circular velocities of the baryons (disk plus bulge) 
and the total mass distribution (baryons plus dark matter) re- 
spectively, both measured at 2.2 disk scale lengths. The solid, 
horizontal line indicates /Scrit: all galaxies above this line con- 
tain a bulge component (whose contribution to V2.2 is given 
by the magenta triangles). Since we tuned (3cnt to match the 
mean bulge-to-disk ratio of the data, its value is different for 
different models. 

In all 3 models, Vgai/V2.2 is strongly correlated with sur- 
face brightness, and to a lesser extent luminosity (middle pan- 
els; see also FAOO and Zavala et al. 2003). In fact, most 
of the luminosity dependence simply owes to the fact that 
more luminous galaxies typically have a higher surface bright- 
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Fig. 10. — Residual correlations for our three models. The upper panels show the residuals of the /-band VL and RL relations, the middle panels show the 
residuals of the stellar mass VM and RM relations, while the lower panels show the residuals of the baryonic mass VM and RM relations. The solid red lines 
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the coiTelation between the residuals decreases in magnitude going from the baryonic mass to stellar mass to 1-band luminosity. 



ness, as is evident from the color-coding. Note that the mean 
Vgai/V2,2 is significantly different for different models. In the 
case of Model 1 (top-heavy IMF), there are no galaxies with 
Vgai/V2.2 > 0.85 (indicated by the gray shaded region). This 
means that, in this model, the baryons never contribute more 
than ^ 70 percent of the total enclosed mass within I.IR^. In 
models II and III the fraction of galaxies with Vgai/V2.2 > 0.85 
are ~ 4% and ~ 37%, respectively. Note that Models I and III, 
though, have a very similar residual correlation slope 7 (up- 
per panels of Fig. [TOb . This is in contradiction to CR99 who 
argued that 7 is uniquely correlated with the mean Vd/V2.2; in 
particular, from the weak residual correlation observed, they 
concluded that high surface brightness disk galaxies should 



have, on average Vd/V2.2 — 0.6. In Appendix ICl we discuss 
the CR99 method in more detail, and show why the observed 
value of 7 does not accurately constrain Vd/V2.2- 

Finally, the lower panels of Fig. [13] show the ratio V2.2/Kir 
as a function of luminosity. Note that the circular velocity at 
2.2 disk scale lengths can be as large as two times the virial 
velocity. The median values of V2.2/Vvir are 1.68, 1.86 and 
1.67 for Models I, II, and III, respectively. As we discuss in 
more detail in 35.1! below, these relatively high values have 
important implications. 
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As discussed the zero point of the RL relation is highly de- 
generate between changes in Agai and mgai.o- In ^we chose 
to tune mgai in order to match the RL relation, while keeping 
Agai fixed at the median value of the spin parameter of dark 
matter halos; Agai = A = 0.042. However, there are several 
reasons why one might expect that the spin parameter of the 
galaxy is different from that of its dark matter halo. Firstly, 
numerical simulations suggest that haloes that are more likely 
to host disk galaxies, i.e., those that did not experience any 
recent major mergers, have a median spin parameter that is 
significantly lower than that of the full distribution of haloes 
(D'Onghia & Burkert 2004; Maccio et al. 2006). Secondly, 
mass is generally more centrally concentrated than specific 
angular momentum. Consequently, if disks form inside out, 
one expects that jgai < mgai, and thus Agai < A. Thirdly, during 
galaxy formation angular momentum may be transferred from 
the baryons to the halo via dynamical friction, also resulting in 
Agai < A. Finally, numerical simulations have shown that the 
directions of the angular momentum vectors of the (hot and 
cold) baryons and the dark matter can be strongly misaligned 
(van den Bosch 2002; Chen et al. 2003; Sharma & Steinmetz 
2005; Bailin et al. 2005b), suggesting that A and Agai are only 



poorly correlated. 

To fully explore these possibilities, we relax the constraint 
that Agai = A = 0.042 and investigate which combinations of 
Agai and mgai match the zero points of the VL and RL rela- 
tions. Our strategy is to run a grid of models in mgai- Agai space 
while keeping all other parameters fixed. In particular, we set 
Q= 1.5, Aimf = -0.2, and /?„!! = 1- For a given combination 
of the concentration normalization, rjc, and the adiabatic con- 
traction parameter, u, we compute the zero point offsets from 
the observed VL and RL relations. We set the zero point to 
logL/ = 10.3. Rather than compute a full model for each Agai 
and njgai, we find the values of Myk that result in logL/ = 10.3. 
This requires an iterative procedure since the stellar mass frac- 
tion, and hence the luminosity, depend on Agai- In order to fa- 
cilitate the interpretation of the zero point offsets we normal- 
ize them to the observed scatter in the VL and RL relations, 
respectively. The results of these grid searches are shown in 
Figim Each row corresponds to a model with a different ry^ 
and u, as indicated. 

Panels in the first and second columns from the left show 
the constraints on Agai and mgai from the VL zero points and 
the RL zero point, respectively. For the RL relation, there is 
a strong degeneracy in all models with Agai oc m'■!;^^ ; however. 
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Fig. 13. — Contribution of the baryons (stellar disk, gaseous disk, and bulge) to circular velocity at 2.2 disk scale lengths, V5ai/V2.2, and ratio of observed to 
virial circular velocity for Models 1, 11, & 111. Point types and colors are the same as in Fig.[T] All models show a strong dependence of Vg^\/V22 with surface 
brightness (upper panels), and a weaker one with luminosity (middle panels). The dependence on luminosity is mostly because more luminous galaxies have 
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contribution of the disk to Vgai is given by the horizontal black solid line. The contribution of the bulge is shown by the magenta triangles, and in all cases 
contributes less than ~ 25% of the total mass within 2.2 disk scale lengths. The lower panels show Vi i/^vk against luminosity. The median V2 2/Kii f™ these 
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for a given mgai the constraint on Agai is very strong. The third 
column from the left shows the quadratic sum of the VL and 
RL offsets, with contours as in the first two columns. The solid 
black dot shows the minimum of the VL + RL offset for each 
model within our grid, so it is not necessarily a global mini- 
mum. The fourth column from the left shows the ratio of the 
circular velocity of the baryons to the total circular velocity at 
2.2 disk scale lengths, Vgai/Va.a- Lines of constant Vsa\/V2.2 are 
approximately Ag^i oc mg^]. In these models Vgai/Vi.a is equiv- 
alent to central surface brightness, so that galaxies with low 
surface brightness disks have low Vgai/Vaa, and vice versa for 
high surface brightness disks. The fifth column from the left 
shows the ratio of V2.2 to Vvh . The dependence of this ratio on 
mgai and Agai is more complicated than Vgai/V2.2 but follows 
the same trend that models with higher mgai and lower Agai 
have higher Vaa/Kh- Thus for model galaxies that match the 
RL zero point, lower Agai solutions have lower Vgai/V2.2 and 
Va.i/Kir- We discuss the significance of this result below. 

The first row uses the standard halo concentrations {rjc = 
1.0) for a ACDM concordance cosmology (On, = 0.3, Ha = 
0.7, ag = 0.9) and standard adiabatic contraction {1/ = 1.0). 
Note that the standard model advocated in MMW, which has 
^Wgai = Agai = 0.05, predicts a VL zero point that is offset from 
the data by ^ 2(t. Matching the VL zero point for these values 
of rjc and ly requires m^ai — Agai — 0.1. However, this clearly 
results in disks that are too large. This demonstrates that this 
'standard' model is unable to simultaneously match the VL 
and RL zero points for realistic T/. 

The second row shows a model with 77^ = 0.5, i.e. with halo 
concentrations that are 50% lower than for the ACDM con- 
cordance cosmology. Although one can now match the VL 
zero point with smaller, more realistic Agai, a simultaneous 
match to the VL and RL zero points still requires relatively 
high mgai. This is essentially Model II. As shown in the third 
row of FiglT4] very similar results are obtained for models 
with standard halo concentrations (rjc = 1) but without adia- 
batic contraction (i^ = 0). Model III is basically an example of 
such a model, which indeed yields results that are very similar 
to those of Model II. 

In the fourth row we consider an even more unorthodox 
model; in addition to 'turning off adiabatic contraction {i/ = 
0), we also set the halo concentrations to be 50% lower than 
expected in the ACDM concordance cosmology (ry^ = 0.5). 
For these parameters, the VL zero point can be matched for 
a very wide range of Agai and mgai. In fact, there are two 
branches in mgai -Agai space that match the VL zero point, cor- 
responding to high and low surface brightness disks. This 
branching is also visible from the curvature of the Agai and 
mgai lines in the upper right panels of Fig. |7] In order to also 
match the RL relation the model needs to fall on the high sur- 
face brightness branch, and needs to have values for both Agai 
and mgai that are much lower than for the previous models. 
Similar constraints are obtained for a model with halo expan- 
sion (v = -l) but with standard halo concentrations (77^ =1), 
shown in the fifth row. 

The lower two rows of panels thus indicate that there is a 
part of parameter space that can match the VL and RL zero 
points equally well as models I-III, but with a much lower 
average spin parameter and much lower galaxy mass frac- 
tions. Before we address the physical relevance of these un- 
conventional models, we investigate to what extent they can 
match the VL and RL slopes and their residual correlations. 
Fig. [15] shows scaling relations and residual correlations for 
one particular model (hereafter Model IV) with halo expan- 



sion and standard concentrations, which matches both the VL 
and RL zero points. The parameters of this model are listed 
in Table 2, for the scatter in the model parameters we adopt 
fine = cTinT = 0.23 (as for models_I-III) and a\n\ = 0.28. The 
median galaxy spin parameter is Agai = 0.023, which is similar 
to the median spin parameter of dark matter halos that have 
not had a recent major merger (D'Onghia & Burkert 2004), 
while mgai.o = 0.03. Both of these values are much lower than 
for models I-III. Note that this model fits the slopes, zero 
points, scatter and correlation between residuals of the VL 
and RL relations. It also predicts a significant correlation be- 
tween the residuals of the corresponding baryonic scaling re- 
lations, in agreement with Models I-III. The average Vgai/Vi.a 
of this model is even larger than for Model III (~ 62% percent 
of all galaxies have Vgai/Vi.a > 0.85), indicating that a larger 
fraction of galaxies is baryon dominated within 2.2 disk scale 
lengths. 

Finally, the lower-left panel of Fig. [15] plots V2.2/Vvir as 
function of the galaxy /-band luminosity. This model pre- 
dicts an average Vaa/Kir of ^ 1 .2, significantly lower than for 
models I-III. Thus, V2.2 in Model IV is comparable to Vmax, the 
maximum circular velocity of a NFW halo without adiabatic 
contraction (see Fig. HI. This is also evident from the panels in 
the fifth column of Fig[T4l which shows contours of constant 
V2.2/Vvir: note that the models in the upper three rows that si- 
multaneously match the VL and RL zero points, indicated by 
a solid dot, all predict relatively high values of Va.a/^vir with 
a mean of ~ 1 .8 ± 0. 1 . However, the models in the lower two 
rows, which can simultaneously fit the VL and RL zero points 
with low Agai and mgai, predict much lower values of V2.2/Vvii- 
with a mean of ~ 1 .2. 

5.1. A revised model for disk formation 

We have shown that one can construct models that match 
the slopes, zero points and residual correlations of the VL 
and RL relations, and which predict relatively low values for 
V2.2/Vvir- This has an important implication: all semi-analytic 
models that are able to match the luminosity function (LF) 
require Vvu < V2.2 ^ Vmax in order to simultaneously match 
the VL zero point (e.g., Somerville & Primack 1999; Cole 
et al. 2000; Benson et al. 2003; Croton et al. 2006). A similar 
conclusion was obtained by Yang et al. (2003) and van den 
Bosch et al. (2003) using the conditional luminosity function 
formalism. This suggests that models I-III will be unable to 
simultaneously fit the observed LF of galaxies. Model IV, 
however, which is based on halo expansion, predicts ratios of 
V2.2/Kir that clearly allow for a simultaneous fit to the LF. 
In more general terms, as can be seen from the fifth row of 
panels in Fig[T4l low values of V2.2/Vvir require average val- 
ues for Agai and mgai that are much lower than the standard 
values normally adopted. In the standard model (i.e., v =\, 
rjc =1), however, such low values for Agai and mgai result in a 
VL zero-point that is in violent disagreement with the data. 
The only models for which a low V2.2/Vvir simultaneously 
allows a fit to both the VL and RL zero points are those in 
the lower two rows: either adiabatic contraction does not oc- 
cur {v = 0) and halo concentrations are a factor two smaller 
than predicted for a ACDM concordance cosmology, or ha- 
los actually expand during disk formation (with :/ ~ -1), in 
which case their original (before expansion) concentrations 
can be as predicted. Of course, one can also construct models 
with intermediate values of v and ijc- For a cosmology with 
(Qm = 0.25, (78 = 0.75) as advocated by the third year WMAP 
results (Spergel et al. 2006) rjc = 0.75 which requires i^ ~ -0.5. 
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Fig . 15. — Scaling relations and residual coirelations of model IV. This model reproduces the slopes, scatter and zero points of the VL and RL relations as well 
as the weak residual correlation. As with models 1, 11, & 111, model IV predicts that the residuals of the baryonic mass VM and RM relations should be strongly 
anti-con'elated. This model contains both baryon and dark matter dominated galaxies at 2.2 disk scale lengths, and has a median V2.2/Kii = 1-22 which is often 
required in order to simultaneously fit the luminosity function and the VL relation. 



This obviously raises the question whether halo expansion 
during disk formation is a realistic option from a physical 
point of view. It is certainly inconsistent with the standard 
picture, in which disks form out of cooling flows that preserve 
their specific angular momentum (Fall & Efstathiou 1980; 
MMW; Dalcanton, Spergel & Summers 1997). Since the time 
scale for new gas to reach the disk is never shorter than the 
free fall time, the formation of the resulting disk will cause 
a contraction of the dark matter halo that is close to adia- 
batic. In order to avoid contraction, one needs to abandon 
the idea that gas reaches the disk via a (relatively smooth) 
cooUng flow. Rather, the disk material has to assemble out 
of several (massive) clumps, which reach the center of the 
halo by dynamical friction. After all, dynamical friction trans- 
fers the potential energy and orbital angular momentum of the 
clumps to the dark matter particles, which consequently move 
to larger radii (e.g., El-Zant, Shlosman & Hoffman 2001; Ma 
& Boylan-Kolchin 2004; Mo & Mao 2004; Tonini, Lapi & 
Salucci 2006). The same mechanism has also been proposed 
as an explanation for the large, constant density cores ob- 
served in some clusters of galaxies (El-Zant et al. 2004; Zap- 
pacosta et al. 2006). 

In addition to providing a scenario for disk formation that 
is consistent with the galaxy luminosity function, halo expan- 
sion may also help to alleviate a long-standing problem with 
observed rotation curves, which are claimed to be inconsistent 
with steeply cusped NEW profiles (e.g., Dutton et al. 2005 
and references therein) and with the observed pattern speeds 
of bars (e.g., Debattista & Sellwood 2000, Weiner et al. 2001). 
Although we defer a detailed discussion on the rotation curve 
shapes of our various models to a forthcoming paper, we have 
verified that in all models discussed here the galaxies have 
realistic (i.e., flat) rotation curves. 



Additional 'support' for Model IV comes from various es- 
timates of the baryonic mass fractions of disk galaxies. As is 
evident from Eig[T4l models with halo expansion that match 
the YL and RL zero points predict much lower values for 
Agai and mgai than the standard models with adiabatic con- 
traction. As we have argiied above, there are several rea- 
sons for expecting Agai < A. Support for low galaxy mass 
fractions comes from various sources, including the condi- 
tional luminosity function formalism (Yang et al. 2003, 2005; 
Cooray 2005) and galaxy-galaxy lensing (Guzik & Seljak 
2002; Hoekstra et al. 2005). In particular, in one of the 
largest galaxy-galaxy lensing studies to date, Mandelbaum 
et al. (2006) find that late-type galaxies with a stellar mass 
of a few times 1O'"/!"'M0, have on average, mgai ~ 0.03 (i.e., 
corresponding to the parameter r\ in Mandelbaum et al. be- 
ing 0.18), in good agreement with the value of mg^i o for our 
Model IV. 

The main potential problem for the 'expansion-scenario' 
proposed here is that it is not clear whether the merger be- 
tween various clumps can produce a realistic disk galaxy. Al- 
though high resolution hydrodynamical simulations of disk 
formation in a cosmological framework have recently sug- 
gested that realistic disk galaxies may form out of merging 
progenitors (e.g., Robertson et al. 2006), more work is clearly 
needed to investigate such a formation scenario in detail. 

To summarize, if one is willing to abandon the idea that disk 
formation involves adiabatic contraction of the correspond- 
ing dark matter halo, one can construct disk formation mod- 
els that simultaneously match the YL and RL relations using 
standard cosmological parameters. In addition, these model 
predict low baryonic mass fractions, in better agreement with 
the data, predict low spin parameters, consistent with a pic- 
ture in which disks form preferentially in halos with quies- 
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cent merger histories, and, most importantly, yield character- 
istic values for V2.2/Vvir that suggest that one may be able to 
simultaneously fit the galaxy luminosity function. 

6. SUMMARY 

We have used the slopes, zero-points, and residuals of the 
observed VL and RL relations to place constraints on the stan- 
dard ACDM -based model for the formation of disk galaxies. 
Our models consider exponential (baryonic) disks in centrifu- 
gal equilibrium within NFW dark matter halos. We model the 
reaction of the dark halo to disk formation in a way which 
permits it to range from the standard adiabatic contraction 
to an expansion of a similar magnitude. The disk stars-to- 
gas ratio is determined by a threshold surface density for star 
formation. A bulge is included based on a self-regulating 
mechanism that ensures disk stability. The disk properties are 
converted into observables using an empirically determined, 
luminosity-dependent, stellar mass-to-light ratio. 

We construct samples of model galaxies including intrinsic 
scatter in halo concentrations, stellar mass-to-light ratios, and 
galaxy spin parameter. In addition, we mimic observational 
errors, and sample the model galaxies so that they reproduce 
the observed luminosity distribution. For each sample we 
construct model VL and RL relations which we compare to 
the data. By demanding the models to simultaneously repro- 
duce the slopes, zero points and scatter of the observed VL 
and RL relations, including the correlation between the resid- 
uals, we obtain the following conclusions: 

• Since the stellar mass-to-light ratio increases with in- 
creasing luminosity, the slopes of the VL and RL rela- 
tions deviate from the basic halo virial relation Vyt oc 

1/3 

Rvir OC M^jj. . The observed VL slope is reproduced 
when the gas-to-stellar mass ratio is properly decreas- 
ing with mass. This is naturally achieved by a surface- 
density threshold for star formation. The observed RL 
slope requires that the disk mass fraction is properly in- 
creasing with mass, in agreement with generic predic- 
tions of supernova feedback (e.g., Dekel & Silk 1986; 
Dekel & Woo 2003). 

• The standard model, assuming adiabatic contraction 
and ACDM halo concentrations, can match the VL zero 
point only with an unrealistically top-heavy IMF. Even 
in this case, the predicted ratio of Va.a/Kii is too high 
for a match of the galaxy luminosity function. 

• Models with a realistic IMF and adiabatic contraction 
require halo concentrations that are a factor two smaller 
than predicted, which is marginally consistent with the 
WMAP constraints on the cosmological parameters. 
Models with realistic IMFs and standard halo concen- 
trations can simultaneously match the VL and RL zero 
points only if adiabatic contraction does not occur In 
both of these cases the predicted ratio of Vaa/Kii- is too 
high for matching the luminosity function. 

• If disk formation causes dark-matter halos to expand 
rather than contract, the VL and RL relations can be re- 
produced with standard halo concentrations combined 
with low values of galaxy spin parameter and mass 
fraction. The low average spin parameter is consis- 
tent with the picture in which disks survive in halos that 
have not experienced recent major mergers (D'Onghia 



& Burkert 2004). The low baryon fraction is consis- 
tent with galaxy-galaxy lensing studies (Mandelbaum 
et al. 2006). This model predicts that V2.2 — 1.2yvi,, 
consistent with what is required to fit the luminosity 
function. 

• The scatter in the VL relation has roughly equal contri- 
butions from scatter in c, in T/ and from observational 
errors. To match the amount of scatter observed, we 
require that a\nc — 0.23. This is smaller than the pre- 
diction for the full set of dark-matter halos, crinc ~ 0.32 
(Bullock et al. 2001a), but consistent with the sub- 
sample of halos without recent major mergers (Wech- 
sler et al. 2002). 

• The RL scatter is dominated by scatter in the spin pa- 
rameter Agai. The observed RL scatter requires that 
CTin A — 0.25, about half the value predicted for ACDM 
halos. This is again consistent with the picture in which 
disks survive in halos with a quiet merger history. This 
picture also implies that the average spin parameter is 
low. A simultaneous match of the VL and RL relations 
with such a low Agai requires that the baryonic mass 
fraction is also low. This, in turn, strongly favors a 
model with expansion rather than adiabatic contraction. 

• The observed residuals of the VL and RL relations show 
only a weak anti-correlation. We have shown that this 
can partly be attributed to the threshold density for star 
formation, which causes the spread in Agai to scatter 
galaxies along the VL relation (see also FAOO). Repro- 
ducing the shallow slope 7 of the residual correlation is 
possible only when there is scatter in T/ and when the 
scatter in c is relatively low, crinc ^ 0.23. 

• Unlike the earlier suggestion by Courteau & Rix 
(1999), we find that 7 does not provide a strong con- 
straint on the baryon fraction in the inner halo. The rela- 
tion between 7 and Vgai/Vai depends on several model 
ingredients such as the degree of halo contraction or 
expansion, the threshold density for star formation, and 
the sensitivity of 7 to the scatter in the different vari- 
ables. In addition, in all our models there is a large scat- 
ter in Vgai/V2.2 which is strongly correlated with surface 
brightness. This correlation results in the highest sur- 
face brightness galaxies being baryon within 2.2 disk 
scale lengths, while at the same time allowing the low- 
est surface brightness galaxies to be dark matter domi- 
nated as is generally accepted. 

• Although observations have shown that the VL and RL 
residuals are uncorrelated, our models robustly predict 
that the residuals of the baryonic VM and RM relations 
are strongly anti-correlated. 

Based on these results we conclude that the standard model, 
which includes standard adiabatic contraction, standard halo 
concentration parameters and a standard IMF, does not allow 
a simultaneous match to the VL and RL relations. This is a 
modification of earlier conclusions, e.g., by MMW and Piza- 
gno et al. (2005), which owes to our more realistic modeling 
(see Appendix IaIi. Although a proper fit can be obtained if 
the halo concentrations are significantly lower, or if adiabatic 
contraction does not occur, these models predict high values 
of V2.2/Vvii which prevent a match with the galaxy luminosity 
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function. They also predict high baryonic mass fractions in 
conflict with results from galaxy-galaxy lensing. 

To circumvent these problems, we advocate a model in 
which the dark matter halo responds by expansion rather than 
by contraction, to the formation of the disk. This model (i) 
predicts values for V2.2/Kir that allow a simultaneous fit to the 
LF (ii) is consistent with disks forming predominantly in dark 
matter halos that have not experienced a recent major merger, 
(iii) predicts low baryonic mass fractions, in agreement with 
a wide range of observations, and (iv) predicts dark matter 
halos that are less centrally concentrated, in better agreement 
with rotation curve shapes and bar pattern speeds. 

The idea of an expanding dark matter halo is counter to the 
standard model for disk formation, in which disks form out 
of relatively smooth cooling flows that conserve their specific 
angular momentum. Rather, the expansion scenario requires 
that disks form out of merging clumps, which transfer energy 



and angular momentum to the dark matter via dynamical fric- 
tion and three body interactions. Such a scenario naturally 
arises if disks form out of clumpy, cold accretion flows (Birn- 
boim & Dekel 2003; Keres et al. 2005; Dekel & Birnboim 
2006). While the formation of realistic disk galaxies out of 
such clumpy streams is yet to be investigated in detail, pre- 
liminary hints are provided by hydrodynamical simulations 
which have demonstrated that disks can originate from merg- 
ing gaseous progenitors (e.g., Robertson et al. 2006). 
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APPENDIX 

A: COMPARISON WITH PIZAGNO ET AL. 2005 

Pizagno et al. (2005, hereafter P05) study the correlations among stellar mass, M,, /-band disk scale length, Rj, and rotation 
velocity at 2.2 disk scale lengths, V2.2, for a sample of 81 disk dominated galaxies (defined to have disk-to-total luminosity 
fractions greater than 0.9). In what follows we refer to the V2.2 -A^* and Ri-M^ relations as the VM and RM relations, respectiv ely. 
P05 estimate stellar masses using g-r colors corrected for extinction and the relations in Bell et al. (2003a). As shown in ^3.41 
the T/-L/ relation of P05 is in excellent agreement with ours (assuming the same IMF). 

P05 claim that a MMW model with m^ = 0.05, A = 0.06, C20() = 10 and using the adiabatic contraction formula of Gnedin 
et al. (2004) matches the VM and RM data reasonably well. This is in disagreement with our statement that the MMW model 
is unable to reproduce the slopes of the VL and RL relations ( ^4. 1.1b and with our claim that models with adiabatic contraction, 
standard halo concentrations, and standard IMFs require ;«gai ~ Agai ~ 0.1 in order to match the VL zero point (Fig.[T4]i. Below 
we address these and other differences between our results and those of P05. 

• We first note that unlike our VL relation and those of Giovanelli et al. 1997, the VM and VL relations of P05 do not follow 
a single power-law. There is a significant deviation below V2.2 = 120kms"' in both their VM and VL relations. In what 
follows we focus on their galaxies with V2.2 > 120kms"'. 

• P05 assume C200 = 10 and that the disk is 100% stars, thus their VM relation has a slope of 1/3, in agreement with their 
data. However, as shown in 34.1.11 the expected variation of c and the gas-to-stellar mass ratio with Mvii both result in 
significantly shallower VM slopes. Thus either c and the gas-to-stellar mass ratio do not depend significantly on Myir, or 
their data significantly over-estimates the VM slope. 

• As can clearly be seen in Fig. 4 of P05, their model with mj = 0.05 and A = 0.06 does not match the VM zero point of the 
data. Their model with mj = 0.10 and A = 0.08 provides a better match, but predicts a much stronger dependence of VM 
scatter on A, and hence disk size. Their models assume pure stellar disks. Including reasonable gas-to-stellar mass ratios 
requires even larger values of m^ and A to match the VM zero point which are now consistent with our results (see Fig. O. 

• P05 do not construct a self consistent model for the VM and RM relations and their scatter. We have verified that pure 
disk models (e.g. MMW) which simultaneously match the VM and RM relations always predict a significant correlation 
between VM and RM residuals. P05 speculated that a weak correlation between size and VM residual could be washed 
out by other sources that contribute to the VM scatter. However the only significant sources to the VM scatter are c and 
observational errors. We have found that while observational errors do reduce the strength of the residual correlation the 
effect is small. By contrast scatter in c results in a significantly stronger negative correlation between the residuals. Thus 
rather than washing out the correlation between size and VM residual, the expected scatter in c and observational errors 
will result in a stronger correlation. 

• The RL relation of P05, as shown in Fig. |3] has a significantly higher zero point normalization than ours. The most likely 
sources of this difference are their bulge-to-disk ratio selection criteria, or their lack of inclination corrections to the disk 
sizes. Their larger disk sizes require a ~ 50% larger spin parameter, Agai, assuming all the other model parameters are kept 
fixed. Larger disks, at a given disk mass, contribute less to V2.2, and result in less halo contraction. This makes it easier to 
fit the VL zero point and also weakens the correlation between the VL and RL residuals. However, as can be seen in Fig. [14] 
even a 50% increase in the RL zero point (corresponding to about Icr of the observational scatter) would not significantly 
change any of our conclusions regarding a simultaneous match the of the VL and RL zero points and a low ratio of V2.2/Vvir- 

In summary, although there are differences in the data samples used by P05 and ourselves, both our and the P05 data require 
large mgai ~ Agai — 0.1 in order to match the VL zero point. Not only is such a high value of Agai unrealistic, such a high njgai is 
also inconsistent with galaxy-galaxy lensing constraints. Furthermore the MMW model adopted by P05 is unable to reproduce 
the slope or the surface brightness independence of the VL relation, due to the overly simplistic assumption that galaxy disks are 
100% stars. 

B : CORRELATIONS BETWEEN MODEL PARAMETERS 

We have assumed that the scatter in the model parameters are uncorrelated. For completeness, we here present a brief discussion 
on how correlations between model parameters might affect our results. 

• A-c anti-correlation: Bullock et al. (2001b) found no correlation between A and c above the weak correlation expected from 
the definition of A. However, Bailin et al. (2005a) claim a correlation, but this is likely due to their inclusion of unrelaxed 
halos (Maccio et al. 2006). Suppose halos with larger A have lower c; looking at Fig. |7]we see that this would result in 
larger scatter in both the RL and VL relations, and by the same reasoning a stronger residual correlation. 

• T/ - c correlation: At a fixed halo mass, lower c halos assemble later; in addition lower c halos result in lower surface 
density disks. Thus, everything else being equal, we expect that at a fixed halo mass, lower c halos to contain galaxies 
with younger (bluer) stellar populations than halos with larger c. Since younger populations correspond to lower stellar 
mass-to-light ratios, we thus expect a positive correlation between c and T/. From Fig. |7]we see that such a correlation 
increases the overall VL scatter It will also decrease the RL scatter, but since T/ and c contribute relatively little to the RL 
scatter, this reduction is unlikely to be significant. 



A Revised Model for Disk Galaxy Formation 27 

• Agai - »Jgai correlation: since mass is more centrally concentrated than specific angular momentum one generally expects 
that jgai < Wgai if disk galaxies form inside out. This implies that a smaller mgai implies a smaller Agai. If these two 
parameters are indeed positively correlated, the overall scatter in the RL relation is predicted to be smaller, as can be seen 
from Fig.|7] However, recall (§4.1 & 4.3) that it is the correlation between the gas-to-stellar mass ratio and Agai, through 
a critical star formation threshold density, that is essential in reproducing the slope of the VL relation and the surface 
brightness independence of the VL relation. Any models that introduce a correlation between Agai and mgai will likely erase 
these successes. 

• T- Agai anti-correlation: At a fixed baryonic mass lower Agai results in higher surface density disks. Since empirically the 
star formation efficiency is proportional to the surface density of the stars (Kennicutt 1998), we thus expect higher surface 
density systems to have older stellar populations, in agreement with the data (Kauffmann et al. 2003b). Such an anti- 
correlation would tend to remove any color dependence of the RL relation and hence reduce its scatter (Bell et al. 2003b). 
However, it would also increase the surface brightness dependence of the VL relation. 

In summary, there are several plausible correlations between model parameters. Although three of these have the potential 
to explain why the RL scatter is smaller than predicted, we do not expect that these correlations can reconcile the full amount 
of scatter in A expected from simulations with the observed scatter in the RL relation. Furthermore, each of these correlations 
will increase the strength of the correlation between VL and RL residuals (or equivalently the strength of the surface brightness 
dependence of the VL relation). Thus to reproduce the weakly correlated residuals we would require even less scatter in c than 
the CTinc = 0.23 we currently adopt. This reinforces the conclusion that disk galaxies form in a subset of halos. 

C: COURTEAU & RIX (1999) REVISITED 

CR99 argued that the weak correlation between the residuals of the VL and RL relations, 7 ~ -0. 1 , implies disks of high surface 
brightness galaxies should have, on average Vd/V2.2 — 0.6. CR99 arrive at this result by computing 91ogV2.2/91og/?d for models 
consisting of exponential disks in NFW halos with adiabatic contraction. They found that in these models dlogVo.i/ d\ogRi 
correlates with Vd/V2.2 (the ratio of the circular velocity of the disk to the total circular velocity, at 2.2 disk scale lengths), such 
that galaxies with higher VijV^.i have more negative t)log V2.2/91og/?d, with the limiting case of 91ogV2.2/91og/?d = -0.5 for a 
pure exponential disk. Thus by assuming that 91ogV2 2/91og/?(i is equivalent to the observed slope of the residual correlation 
7, they concluded that yd/^2.2 — 0.6. This would suggest more DM than baryons within 2.2 disk scale lengths, but not that the 
baryon contribution is insignificant. If the baryons were insignificant then the CR99 model would predict 91og V2.2/91og/?d > 0. 
However, as outlined below, several of the assumptions made in CR99 are not true in general. 

• CR99 assumed that halos contract adiabatically to the formation of the disk. In Fig. [16] which should be compared with 
Fig. 9 of CR99, we repeat the CR99 analysis for halos with and without adiabatic contraction. Note that CR99 used 
C200 = 10 and that Cvir — 1.3c2oo, for ease of comparison we also adopt their definition here. We see that turning off 
adiabatic contraction increases Vd/V2.2 by about 30% at a fixed c)logV2.2/91og/?d- 

• CR99 assumed that bulge formation does not affect the scale length of the disk. We have shown that in models where 
the specific angular momentum of the disk increases due to bulge formation (e.g. /x = 0.25) the residual correlation is 
weakened. 

• CR99 assumed that the contribution of the gas to the disk is negligible. While this is true for t he highest surface brightness 
galaxies, moderate surface brightness galaxies contain more gas. We have shown (see ^4.3b that the correlation between 
the gas-to-stellar mass ratio and disk surface density significantly reduces the strength of the residual correlations. 

• CR99 underplayed the effect of scatter in halo concentration c and stellar mass-to-light ratio T. Fig.[T2]shows that scatter 
in these parameters can result in significantly different residual correlation slopes. Scatter in c results in a more negative 
correlation while scatter in T/ results in a more positive correlation. For realistic amounts of scatter in both c and T/ the 
residual correlation is weakened. 

• CR99 assumed that 7 = 91ogV2.2/91og/?d- In general this is not true, because 7 is a global quantity while 91ogV2.2/91og/?d 
is a local quantity. More specifically 7 is the slope of the VL and RL residuals for a sample of galaxies, while 
(91ogV2.2/c'log7?d is the slope of the change in V due to change in R, at fixed disk mass, for a single galaxy. 

Thus we conclude that the VL-RL residuals cannot be used to place model independent constraints on the baryonic fraction of 
disk galaxies. However, this does not diminish their importance as a constraint for galaxy formation models. 
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Fig. 16. — CR99 analysis for NFW halos with C200 = 10 with and without adiabatic contraction (AC). Note that c,.ir — 1.3c2oo. The y-axis shows the change 
in logV2.2 due to changes in logRj at fixed Mj. The line and color types correspond to different values of V200 and Rj, respectively, as indicated in the figure. 
The horizontal solid and dotted lines show the observed slope of the AlogV- Alogi? relation. The vertical dotted lines show the Vd/V2.2 corresponding to the 
V200 = 180,i?(j = 4 models. Models without AC have a substantially larger Vd/^2.2 (about 30%) than models with AC. 



